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Exact,  closed  form  solutions  for  diffraction  and  backscattering  by  a 
half-plane  whose  surfaces  exhibit  arbitrary  Impedances  were  determined  for 
the  cases  of  a point  source,  a line  source,  and  a plane  wave.  An  Integral 
method  was  devised  that  could  account  for  the  Impedance  boundary  conditions. 
Since  the  solution  is  in  a closed  form,  interpretation  of  the  influece  of  the 
Impedance  covered  half-plane  on  the  diffraction  field  is  greatly  facilitated. 
Specifically,  the  effect  of  each  surface  on  the  acoustic  field  dominates  in  tie 
half-space  in  which  it  faces.  Though  the  effect  of  each  surface  extends  into — ^ 
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Is  very  small  and  diminishes  even  more  as  the  observer  moves 
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CHAPTER  I 


INTRODUCTION 


1 . 1 Background 

Whenever  Che  propagation  of  a sound  wave  is  obstructed  by  some 
obstacle,  part  of  Che  sound  wave  is  deflected  or  disturbed.  This 
disturbance  or  interference  gives  rise  Co  Che  phenomenon  called 
diffraction. 

The  problem  of  diffraction,  where  Che  interfering  obstacle  is 
a semi-infinite  half-plane,  is  a classical  problem  in  Che  area  of 
wave  propagation.  Recently,  because  of  environmental  concerns, 
there  has  been  interest  in  using  barriers  and  enclosures  as  a means 
of  isolating  noise  from  adjacent  areas.  Understanding  the  phenomenon 
of  diffraction  is  important,  since  Che  sound  waves  will  diffract 
around  the  barrier.  Although,  the  geometry  of  Che  half-plane  plays 
a dominant  role  in  creating  Che  diffracted  pressure,  Che  nature  of 
Che  surface  Impedance  also  Influences  Che  behavior  of  the  diffracted 
pressure.  By  adjusting  the  surface  Impedance  of  Che  half-plane, 

Che  diffracted  pressure  could  be  modified  in  some  desirable  fashion. 

1.2  Review  of  Previous  Investigations 

Studies  of  diffraction  have  been  extensively  carried  out  in 
the  fields  of  acoustics  and  optics.  Many  related  publications  may 
be  easily  found  in  the  literature.  This  review  is  concerned  specif- 
ically with  chose  papers  chat  are  related  to  Che  diffraction  of  a 
wave  by  a half-plane  and,  in  particular,  an  Impedance  covered  half- 


plane. 
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Because  of  the  geometry  of  the  half-plane,  the  usual  techniques 
for  solving  boundary  value  problems  do  not  apply.  Consequently,  in 
order  to  solve  the  problem  for  sharp  geometries  like  the  half-plane, 
other  methods  have  had  to  be  formulated.  It  was  Soramerfeld  [1,2]  who 
successfully  used  an  Integral  expression  to  represent  the  scattered 
pressure.  With  the  pressure  expressed  in  an  Integral  form  and  a 
concept  of  sources  and  image  sources,  he  matched  the  conditions  on 
Che  surfaces  of  the  half-plane.  In  Sommerf eld ' s classical  work,  he 
was  able  to  solve  Che  problem  only  for  Che  rather  ideal  boundary 
conditions  of  a perfectly  soft  (pressure  release)  or  a perfectly 
hard  (rigid)  surface.  Any  impedance  between  Chose  extremes  could 
not  be  accommodated  by  his  more  or  less  heuristic  formulation  of  the 
problem.  Sommerfeld's  solution  accounts  for  the  geometry  of  the 
diffracting  obstacle  but  is  limited  in  the  sense  that  it  does  not 
account  for  arbitrary  surface  conditions. 

After  Sommerf eld,  several  researchers  found  a number  of  ways 
to  solve  the  problem  of  diffraction  by  a half-plane  based  on  the  so- 
called  Welner-Hopf  Technique.  The  Weiner-Hopf  method  is  essentially 
an  extension  of  Four ler-Laplace  transforms  which  results  In  a complex 
variable  equation  that  is  analytic  within  a certain  region.  The 
fundamental  step  in  Che  Welner-Hopf  method  is  to  determine  two 
functions,  each  of  which  satisfy  a complex  variable  equation  in 
different  halves  of  the  complex  plane  but  together  satisfy  the 
equation  in  a proper  overlapping  region.  The  problem  can  b?  reduced 
to  the  task  of  splitting  a complex  function  into  two  factors  in  such 
a way  that  each  factor  contains  specific  singularities.  A variation 
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of  the  Weiner-Hopf  method  is  a dual  Integral  technique  employed  by 
Jones  [3]  and  Clemmow  [4].  In  this  method,  two  integral  equations 
are  found  that  are  valid  in  separate  half-planes  of  the  complex 
plane.  It  is  necessary  to  find  a function  that  will  satisfy  each 
Integral  equation  in  its  respective  half-space.  Jones  and  Clemmow 
applied  the  dual  Integral  method  to  the  specific  problem  of  diffrac- 
tion by  a half-plane  and  found  a solution  exactly  equivalent  to 
Sommerfeld' s solution  Again,  they  only  considered  ideal  boundary 
conditions. 

The  Welner-Hopf  technique  was  successfully  employed  by  Senior 
[5].  He  attempts  a problem  in  which  a semi-infinite  metallic  sheet 
has  some  finite  conductivity.  The  solution  that  he  finds  is  exact 
but  the  function  resulting  from  the  Welner-Hopf  analysis  Is  in 
Integral  form  and  cannot  be  easily  computed. 

Following  Senior,  Williams  [6]  extended  the  Weiner-Hopf  technique 
to  a more  general  problem  of  the  diffraction  by  a wedge.  He  reduced 
the  problem  to  the  solution  of  an  ordinary  difference  equation.  This 
difference  equation  was  then  solved  by  means  of  a double  Gamma 
function  defined  by  Barnes  [7].  If  the  wedge  angle  was  of  the  form 
p 7r/q,  where  p and  q are  prime,  and  p is  odd,  a closed  form  solution 
could  then  be  found.  By  choosing  proper  values  of  p and  q,  a wedge 
that  is  very  close  to  a half-plane  could  be  approximated.  Rawlins 
[8]  used  the  Welner-Hopf  technique  to  solve  Che  diffraction  by  a 
half-plane  with  one  side  soft  and  the  other  side  rigid.  According 
to  Rawlins,  a "perfectly  absorbing"  surface  can  be  modeled  by  adding 
the  hard  and  soft  surface  solutions  together  and  dividing  by  two. 
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t From  this  he  concludes  that  the  absorbing  surface  will  attenuate  the 

diffracted  pressure  If  It  Is  on  the  Illuminated  surface. 

I The  case  of  a half-plane  with  an  arbitrary  surface  Impedance 

has  not  been  solved  by  the  Welner-Hopf  technique.  However,  the 
Russian  physicist,  Malyuzhlnets  [9,10,11]  has  approached  the  problem 
by  yet  another  method.  Malyuzhlnets  used  a Fourier  Transform 
approach  to  find  a single  spectral  function  for  the  whole  range  of 
Integration.  The  function  Is  found  to  be  a product  of  four  integral 
functions  that  comprise  the  special  functions  derived  by  Malyuzhlnets 
to  account  for  the  arbitrary  surface  impedance.  Malyuzhlnets 
j demonstrates  that  each  spectral  function  is  a solution  of  a differ- 

ence equation.  Malyuzhlnets  and  his  colleagues  have  also  applied 

I 

^ their  technique  to  a more  general  class  of  problem  where  the  dif- 

fracting obstacle  Is  an  Impedance  covered  wedge. 

In  most  cases,  the  diffraction  problem  was  solved  considering 
an  Incident  plane  wave.  Redfearn  [12]  considers  various  kinds  of 
sources.  Specifically,  for  a point  source  Redfearn  discussed  the 
practical  uses  of  a half  plane  as  a noise  barrier.  Finally,  Clemmow 
[4]  discusses  the  diffraction  of  a line  source  by  a half-plane, 
where  he  derived  the  solution  In  terms  of  the  Fresnel  Integrals.  A 
I good  summary  of  the  work  accomplished  in  Che  area  of  diffraction  by 

I edges  and  wedges  is  given  by  Pierce  [13]. 


CHAPTER  II 


MATHEM/UICAL  FORMULATION  FOR  THE  DIFFRACTION  OF 
ACOUSTIC  WAVES  BY  AN  IMPEDANCE  COVERED  HALF-PLANE 

2 . 1 Statement  of  the  Problem 

In  this  study,  the  problem  of  diffraction,  where  the  disturbing 
obstacle  Is  a half-plane.  Is  considered.  The  acoustic  Impedance  of 
the  upper  and  lower  surfaces  of  the  half-plane  will  have  an  effect 
on  the  diffraction  of  a sound  wave.  Conventionally,  this  problem 
Is  approached  by  considering  only  the  Ideal  surface  conditions  of  a 
perfectly  rigid  or  a perfectly  soft  half-plane.  In  this  study, 
however,  the  effect  of  more  realistic  surface  conditions  are  con- 
sidered. These  surface  conditions  can  account  for  an  arbitrary 
local  reactance  and  an  arbitrary  absorption  on  both  surfaces  of  the 
half-plane.  Specifically,  this  study  alms  to  determine  exactly  how 
the  surface  impedances  on  the  two  sides  of  a half-plane  will  affect 
the  diffraction  of  sound. 

In  addition  to  this,  various  kinds  of  acoustic  sources  can 
produce  the  sound  waves  Incident  upon  Che  half-plane.  Subsequently, 
this  study  also  Investigates  Che  diffraction  for  a plane  wave,  a 
line  source,  and  a point  source. 

Often  when  dealing  with  the  diffraction  of  a sound  wave, 
interest  concentrates  only  on  how  the  diffracted  field  appears  when 
an  observer  Is  located  many  wave  lengths  from  the  diffracting 
obstacle  However,  it  Is  a goal  of  this  study  to  also  investigate 
Che  diffraction  when  Che  observer  distance  from  the  edge  may  not 
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necessarily  be  great,  as  well  as  when  the  observer  distance  does 
become  large. 

A special  case  of  diffraction  results  when  the  source  and 
observer  are  located  at  the  same  position.  This  Is  called  back- 
scattering.  Results  of  backscatterlng  for  various  surface  and  source 
conditions  are  also  to  be  analyzed. 

In  short,  the  goals  of  this  study  are  to  develop  analytical 
models  and  exact  solutions  that  account  for  arbitrary  realistic 
surface  conditions,  source  and  observer  locations,  for  the  acoustic 
diffraction  and  backscatterlng  by  a half-plane. 

2 2 Assumptions 

The  problem  that  Is  modeled  here  Is  based  upon  the  existence 
of  a seml-lnf Inite  thin  plane.  Needless  to  say,  one  cannot  realize 
such  an  edge  physically. 

Perhaps  of  more  significance  are  the  assumptions  and  restric- 
tions that  are  made  concerning  the  acoustical  behavior  of  the 
semi-infinite  plane  surfaces.  For  several  reasons,  the  surfaces 
are  considered  to  have  locally  reacting  impedances.  This  is  chosen 
In  the  first  place  because  it  represents  many  real  acoustic  situations 
reasonably  well.  Secondly,  it  can  be  handled  mathematically  by  a 
comparatively  simple  formulation. 

What  then  Is  meant  by  a locally  reacting  surface?  Suppose  an 
acoustic  wave  impinges  on  a surface.  Some  of  the  energy  of  the 
acoustic  wave  could  be  reflected,  some  could  be  absorbed,  and  some 


could  be  used  to  set  the  surface  Itself  Into  motion.  If  the  surface 
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Itself  Is  set  Into  motion,  It  Is  possible,  depending  on  the  material 
making  up  the  surface,  for  that  motion  to  propagate  along  the  surface 
For  a locally  reacting  surface,  this  Is  not  possible. 

It  may  be  possible  for  the  surface  to  be  set  into  motion  by  an 
Incident  acoustic  wave  In  such  a way  that  each  point  on  the  surface 
moves  independently  of  all  adjacent  points.  In  other  words,  there 
Is  no  coupling  of  the  motion  of  the  neighboring  points  of  the  surface, 
so  consequently,  no  allowances  are  made  for  the  propagation  of  a 
wave  along  or  through  the  surface.  It  Is  Important  to  recognize 
that  this  means  there  are  no  elastic  waves  propagating  In  the  surface 
nor  are  there  any  waves  transmitting  through  the  surface.  However, 
surface  waves  In  the  acoustic  medium  are  accounted  for  In  this 
analysis. 

The  motion  of  the  surface  at  a given  point  Is  determined  by  the 

pressure  of  the  acoustic  wave  at  that  point.  It  turns  out  that  the 

pressure  Is  proportional  to  the  component  of  velocity  normal  to  the 

surface.  The  motion  of  the  acoustic  fluid  tangent  to  the  surface 

will  not  cause  any  motion  in  the  surface.  This  is  the  true  provided 

the  acoustic  fluid  Is  considered  Invlsld  or  at  least  has  negligible 

viscosity.  Thus,  the  ratio  of  the  acoustic  pressure  at  the  surface 

to  the  outward  normal  component  of  the  velocity  at  the  surface  at 

a given  point  Is  called  the  acoustic  impedance  of  Che  surface  at 

that  point.  It  is  usually  signified  by  Z ■ ^ . The  half-plane  is 

n 

considered  to  have  different  impedances  on  Its  two  surfaces  which 
will  be  denoted  by  and  Z , respectively. 
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In  general,  Che  acoustic  impedance  of  Che  surface  may  be  complex. 
The  Imaginary  part  of  Che  Impedance  corresponds  Co  Che  reactive 
component  of  the  surface  Impedance.  This  means  Chat  Che  surface  can 
be  locally  idealized  as  consisting  of  a series  of  isolated  masses 
and  springs.  The  real  part  of  the  Impedance  physically  represents 
the  acoustic  energy  that  is  absorbed  by  Che  surface. 

In  brief,  Che  assumptions  Imposed  to  facilitate  the  mathematical 
formulation  are  (1)  a thin  half-plane,  and  (ii)  locally  reacting 
Impedances  on  Che  surfaces. 

2 . 3 Integral  Representation 

In  order  to  have  an  exact  solution  to  a diffraction  problem,  it 
is  necessary  Co  satisfy  three  basic  requirements.  First,  Che  solu- 
tion must  satisfy  Che  wave  equation.  Second,  the  solution  must 
account  for  the  boundary  conditions  on  Che  surfaces  of  the  diffracting 
half-plane  and  be  continuous  throughout  the  physical  space.  Third, 
the  solution  must  correspond  to  the  given  type  of  excitation,  be  it 
a plane  wave,  a line  source,  or  a point  source.  If  all  these  require- 
ments are  met,  then  Che  resulting  solution  must  be  unique  to  the 
problem  as  it  is  posed.  In  this  regard,  Che  analytical  method  used 
in  this  study  exploits  an  integral  representation  that  must  satisfy 
these  requirements. 

Since  the  Fourier  Integral  transforms  are  useful  in  the  formu- 
lation of  the  problem,  they  are  introduced  here.  Define  the  Fourier 


Transform  pair  as 
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f(x)  - j F(Y)e‘^'*  dY 


(2.1) 


and 


f(y) 


f 


£(x)e^^’‘  dx 


(2.2) 


Consider  now  the  wave  equation  for  the  pressure 


„2  1 a p 

^ P ■-  “2 

C dt 


(2.3) 


where  c Is  the  velocity  of  propagation  of  an  acoustic  wave  In  free 
space.  Assuming  the  pressure  to  be  periodic  In  time  as  e^^^ , then 
Equation  (2.3)  becomes 


V^p  + k^p  - 0 , (2.4) 

where  k - w/c. 

Applying  Fourier  Transforms  on  Equation  (2.4),  one  obtains 


^ + (k^  - Y^)  P - 0 
dy 


(2.5) 


or 


^ + 6^  - 0 

dy 


(2  6) 
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where 


The  solution  to  the  differential  equation  in  Equation  (2.6)  is 

p - A(-y)e'^^^  . (2.7) 


The  (+)  or  (-)  sign  in  the  exponent  corresponds  to  the  two  possible 
independent  solutions  of  Equation  (2.6).  To  choose  the  proper  sign, 
consider  the  physical  situation.  The  pressure  p is  meant  to 
represent  the  diffraction  from  the  edge.  Consequently,  the  solution 
must  reflect  the  behavior  of  outgoing  waves  from  the  edge  rather 
than  an  Incoming  wave  from  infinity.  The  requirement  for  outgoing 
waves  can  be  assured  by  using  the  plus  sign  (+)  i.or  y<0,  and  the 
minus  sign  (-)  for  y>0. 

Taking  the  inverse  transform  of  p gives  the  following  integral 
representation , 

p(x,y)  - ^(Y)e"^®^  e dy  (2.8) 

Since  the  integral  representation  in  Equation  (2.8)  satisfies 
the  wave  equation,  p(x,y)  is  then  chosen  as  a suitable  representation 
for  the  scattered  pressure.  Letting  A “ y/V.  then  Equation  (2.8) 


becomes 
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p(x,y) 


_k 

2tt 


A(A)e 


±iky/l-X^ 


-IkAx 

e 


dA 


(2.9) 


Another  way  Co  consider  this  Integral  representation  Is  through 
a more  heuristic  explanation  of  the  physical  meaning.  The  exponent 
of  the  Integrand  can  be  regarded  as  waves  propagating  in  all  direc- 
tions in  Che  x,y  plane.  This  part  will  satisfy  Che  wave  equation. 
The  function  A(A)  of  the  integrand  can  be  regarded  as  a spectrum 
function.  This  is  to  say  Chat  A(A)  represents  a weighting  or 
distribution  of  waves  that  will  propagate  in  any  given  direction. 
From  this  explanation,  it  is  argued  that,  since  Che  integral  itself 
is  made  up  of  a superposition  of  waves  that  individually  satisfy 
the  wave  equation,  then,  this  superposition  of  waves  in  the  form  of 
the  Integral  will  also  satisfy  the  wave  equation.  Subsequently,  the 
Integral  in  Equation  (2.9)  is  often  referred  to  as  a spectrum  of 
plane  waves. 

A final  step  at  arriving  at  an  appropriate  integral  representa- 
tion is  CO  express  the  integral  of  Equation  (2.9)  in  polar 
coordinates.  This  is  accomplished  by  adopting  two  transformations. 
First,  let 


X - r cos$, 

and 


y “ r sln(J)  . 


(2  10) 
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Next  use  Che  following  change  of  variable  for  X: 

X - cosa  . (2.11) 

Thus,  Che  corresponding  path  of  Integration  In  Che  cx-plane  must 
be  computed. 

Making  use  of  Equations  (2,10)  and  (2.11),  Che  Integral 
expression  may  be  expressed  In  polar  coordinates  as  follows: 


p(r,4)) 


A(cosa)  slna 


^-lkr(cosi(i  cosu  t sln^i  slna)^^ 


(2  12) 


with  (+)  for  the  upper  half-plane  when  y>0,  and  (-)  for  the  lower 
half-plane  when  y^O,  and  where  y Is  the  path  of  Integration.  Letting 
P(cosa)  * k A(cosa)  sina  and  using  the  appropriate  Identity  In  Che 
exponent,  one  arrives  at  Che  following  integral  representation: 


p(r,$) 


P(cosa) 


-ikr  cos((J>+a)  . 

I da 


(2.13) 

(-)  upper  half-plane 
(+)  lower  half-plane 


In  this  form.  It  Is  easier  to  recognize  that  Equation  (2.13) 
represents  an  angular  distribution  of  plane  waves. 

It  Is  necessary  to  carefully  investigate  what  occurs  under  the 
change  of  variables  chosen  In  Equation  (2.11).  Starting  with  the 
equation  for  A: 


X - cos(a^  + 1 a^) 
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and 


A + i A.  - cosa  coslioi^  - 1 sinha^  pina  , (2.14) 

rl  rl  i r 

» 

where  A and  a are  Che  real  parts  of  A and  a and  A.  and  a,  arc  the 
r r 1 i 

Imaginary  parts  of  A and  a.  Equating  the  real  and  imaginary  parts, 
it  follows  Chat 

A “ cosa  cosha^  (2.15) 

r r i 

and 

A^  “ - sinha^  sina^  . (2.16) 

Due  to  Che  periodicity  of  cosa^  and  sina^,  the  entire  complex 
A-plane  will  map  into  a strip  of  width  tt  in  Che  a-plane,  and  Chen 
repeat  itself  in  subsequent  strips  of  width  equal  Co  tt. 

To  understand  this,  consider  A^  negative,  or  Che  lower  half  of 
the  complex  A-plane,  then  - sinhaj  slna^  must  be  negative.  This  is 
met 

1)  for  a^<  0 and  sina^<  0 
if  (2n-l)Tr  £ a £ 2m 

and 

2)  for  a.>  0 and  slna  > 0 

1 r 


if  2nTT  < a < (2n+l)''i 


(2.17) 
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If,  on  the  other  hand,  Is  positive  and  thus  represents  the  upper 
half  of  the  complex  X-plane,  then  - sinhu^  slnu^  must  also  be 
positive.  Consequently,  this  condition  is  also  satisfied: 


1)  for  a < 0 and  sinu  > 0 

i r 

if  2nTT  < a < (2n+l)iT 
— ir- 

and 

2)  for  a^>  0 and  8inoi^<  0 

if  (2n-l)iT  < o < 2nTr  . (2  18) 


Using  the  conditions  stated  in  Equations  (2.17)  and  (2.18), 

Figure  2 1 shows  the  mapping  of  the  entire  complex  X-plane  to  a 

multiplicity  of  strips  In  the  complex  a-plane. 

Because  the  X-plane  maps  Itself  every  n in  the  a-plane,  it  then 

follows  that  the  path  of  integration  from  - “ to  + “ in  the  X-plane 

must  likewise  exhibit  the  same  repetitive  behavior.  Considering  the 

expression  for  the  real  part  of  X in  the  Equation  (2.13)  to  determine 

the  path  for  X from  - <»  to  + “>,  certain  values  for  a and  a,  must  be 

r i 

chosen.  Starting  at  - » and  moving  to  + « in  the  X-plane  the 
corresponding  values  for  and  a^  are  Indicated  in  Table  2,1.  From 
Table  2.1,  it  follows  that  the  path  of  integration  in  the  a-plane 
can  be  made  by  picking  any  combination  of  values  of  a^  and  a^. 
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f 

Table  2,1.  The  path  of  integration  in  the  X-plane  and  the 
Corresponding  a-plane. 


X 

— 00 

-1  -► 

0 - 

1 -► 

OO 

±co 

0 

0 

0 

+O0 

a 

r 

(n+l)Ti 

(n+l)7T 

n+1 

nTT 

nv 

There  Is  one  last  detail  to  consider  before  actually  defining 
\ the  Y-path.  There  Is  a negative  sign  introduced  in  the  integral  by 

the  change  of  variable  In  Equation  (2.11).  To  compensate  for  that 
minus  sign,  Che  path  Is  simply  followed  In  the  opposite  direction. 
Taking  this  fact  into  account,  and  the  values  of  and  from 
Table  2.1,  the  possible  paths  that  may  define  the  ypath  in  the 
a-plane  are  shown  in  Figure  2,2.  Any  path  that  Joins  A - + » to 
X - - 00  is  an  acceptable  Y”path  for  the  integral  representation  of 
Equation  (2.13),  It  is  now  evident  that  there  is  also  a multiplicity 
of  Y-paths. 

Although  all  the  Y~paths  described  in  Che  previous  section  do 
indeed  represent  a path  from  - oo  to  + oo,  not  all  paths  are  suitable 


choices.  The  reason  lies  in  the  behavior  of  the  exponent  in  the 
integrand  of  Equation  (2.13).  One  requirement  for  this  integral  to 
exist  is  that  it  mucL  converge.  Since  the  integration  is  in  the 
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complex  a-plane,  one  must  be  careful  to  Insure  that  the  Y~p3th  will 
lie  In  a region  where  the  exponent  will  insure  convergence  of  the 
integral.  Thus  one  can  find  these  regions  in  the  a-plane  by  finding 
where  the  real  part  of  the  exponent  is  negative  as  follows: 


Ke 


■ikr 


/, 


cos(())-a^)  cosha^  ± i 


slnha.  sln((()- 


“r>}]  ' "■ 


(2.19) 


Consider  the  case  when  y > 0,  so  one  may  choose  the  upper  sign. 
Then  the  real  part  of  the  exponential  becomes  k.rsinha^sln((ti-a^) 

For  a^<  0,  then,  sln((})-a^)  > 0 which  Implies 


$ + (2  -1)71  ^ £ 0 + 2 7T 


For  a^>  0,  then,  sln(4i-a^)  < 0,  so  that 


+ 27T<a<(|i  + (2.  +1)7^ 


For  any  arbitrary  angle  i})  in  the  upper  half-plane,  the  regions  of 
convergence  with  a suitable  v-path  are  shown  in  Figure  2.3.  The 
regions  of  convergence  are  of  width  77  and  also  exhibit  repetitive 
behavior 

The  choice  of  the  Y~pa^b  also  will  be  sufficient  to  satisfy 
the  Sommerf eld  Radiation  Condition.  For  the  two-dimensional  case 
of  plane  wave  radiation,  this  condition  is 
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r 

Substituting  the  expression  for  the  scattered  pressure  into  the 
Sommerfeld  Radiation  Condition  yields 


+ ikp 


(2.20) 


lira 


P ( CO sa)  (1-cos  (iJ)Ta) ) 


-ikr  cos(4)+a) 
e 


0 


(2.21) 


Again,  this  limit  can  be  satisfied  if  the  y-path  lies  in  a region 
where  the  exponential  decay  allows  convergence  as  r approaches 
infinity-  Thus,  the  Sommerfeld  Condition  is  satisfied  by  the  proper 
choice  of  y-path. 

2 . U Formulation  of  the  Integral  Equations  for  an  Incident  Plane  Wave 

The  preceding  section  shows  that  an  integral  representation  for 
the  scattered  pressure  does  comply  with  one  of  the  basic  requirements 
for  a solution.  Specifically,  it  satisfies  the  homogeneous  wave 
equation  and  the  Sommerfeld  Radiation  condition  through  the  choice 
of  y-path.  In  this  section,  the  details  of  applying  this  representa- 
tion to  the  diffraction  of  the  half-plane  with  impedance  boundary 
conditions  will  be  worked  out. 

The  diffracting  obstacle  is  the  half-plane  defined  by  x < o, 
y " o and  - <»  < z < + <».  The  integral  representation  to  be  used  Is 
expressed  in  terms  of  polar  coordinates  so  one  measures  angles  from 
the  positive  x-axls.  Thus,  the  upper  surface  of  the  half-plane  is 
located  at  an  angle  (f)  • + v,  while  the  lower  surface  is  located  at 


A 


an  angle  (()*•-  tt.  An  incident  plane  wave  can  be  represented  as 
follows : 


ikr  cos((t-<b  ) 

P 6 O 

O 


(2.22) 


The  angle  (Ji^,  also  measured  from  the  positive  x-axis , will  denote 
the  angle  of  incidence.  Finally,  the  upper  and  lower  surfaces  have 
impedances  of  and  Z , respectively.  The  geometry  for  the 
half-plane  it  depicted  in  Figure  2.4. 

The  surface  impedance  is  the  ratio  of  the  total  pressure  at 
the  surface  to  the  outward  normal  component  of  the  total  velocity 
at  the  surface.  The  transverse  component  of  velocity  may  be  ex- 
pressed in  terms  of  the  pressure  p by 


V 

0 


. 1 -E 
kpc  r 30 


(2.23) 


It  follows  that  the  transverse  component  of  the  velocity  for  the 
incident  plane  wave  is  given  by 


y _ 1 3 ikrcos(0-0  ) 

i0  kpcr  30  ^o®  °J 


p^sin(0  0^/  ikrcoa((ti-0  ) 

e o 

pc 


(2  24) 


and  also  the  velocity  for  the  diffracted  pressure  is  given  by 


OBSERVER  (r,<^) 


SOURCE (^o) 


Figure  2.4.  Geometry  for  the  impedance  covered  half-plane. 
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(f)  kpcr  3<{) 


i 


P(cosa)  e 


■Ikrcos  (iji+a) 


da 


_sin(c^  P(co8a)  e-lkrcoB(^*cx)  _ 

pc 


(2.25) 


The  following  boundary  conditions  at  the  upper  and  lower  faces  of 
the  half-plane  must  be  satisfied, 


P 

tot 

V 

tot 


Pj+P 

(v^+v) 


±TI 


+ 

z 


(2  26) 


or 


Pj^+P 


±TT  ■ ='"<-1-^'') 


±1T 


(2.27) 


Using  the  expressions  found  in  Equations  (2,13),  (2.22)  (2.23)  and 
(2.24)  for  the  pressures  and  velocities,  and  from  the  boundary 
condition  of  Equation  (2.27),  one  can  determine  specific  relationships 
that  satisfy  the  boundary  conditions.  The  boundary  condition  for 
the  upper  surface  at  4)  ••  tt  becomes: 


, ikrcosa.  . -lkrcos<(in 
P(co8a)e  da+p^e  o 


+ r 

z 

pc  ‘ 


^ Ikrcosa,  , ^o*  -ikrcos0. 

sinaP(co8a)e  da+  : e 0 . 


pc 


(2.28) 
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A similar  expression  for  the  lower  surface  at  ({)  " - n'  is 


, ikrcosa,  ^ -ikrcos4>o 
P(cosa)e  da  + p e 


^ f 

pc  I 


sinaP(cosa)e^‘^'^‘^°«°‘da  - sln^  ^-ikrcos4>o  ^2.29) 

pc  o 


Let 


3ln6 


+ pc 


and 


sine  - , (2.30) 

z 

where  9 and  0 are  known  as  Brewster  angles.  These  angles  have  an 
important  physical  significance.  There  will  be  no  reflection  of  a 
plane  wave  from  an  infinite  plane  if  the  angle  of  incidence  of  the 
plane  wave  is  equal  to  the  Brewster  angle.  In  short,  it  is  the 
angle  of  no  reflection. 

Using  the  relationships  in  Equation  (2.30),  Equations  (2  28) 
and  (2,29)  may  be  rewritten  as  follows: 


/, 


P(cosa)  (8in9^+8lna)e^*^*’*'°®^da  " - p (8ln6^-8ine  )e 

o o 

(2.31) 


d 
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and 


( P (cosa)  I'sinG  -sina)e^*^*'*'°^'^da  ■ - p (sin6  +sin<)  )e 
I o 

(2.32) 

Equations  (2.31)  and  (2  32)  represent  two  Integral  equations  that 
satisfy  the  Impedance  boundary  conditions  - 

It  Is  necessary  that  the  Integral  representations  for  the 
pressure  and  the  velocity  be  continuous  In  the  physical  space  There 
Is  a change  In  sign  In  the  exponent  of  the  Integral  representation 
as  one  crosses  from  the  upper  half-plane  to  the  lower  half-plane. 

To  guarantee  continuity,  the  pressure  and  velocity  In  the  upper 
region,  y > 0,  and  the  pressure  and  velocity  in  the  lower  region, 
y < 0,  must  be  equal  where  the  two  regions  are  joined;  namely,  at 
ifi  “ 0°.  The  two  regions  are  not  joined  together  at  it  or  -tt  because 
the  half-plane  physically  separates  them.  Checking  the  continuity 
of  the  pressure  at  C « 0°,  it  follows  from  Equations  (2.13)  and  (2  22) 
that 


, -ikrcosa.  . ikrcos<tin 
P(cosa)e  du  + p e ° 


. -ikrcosa,  , ikrcoscjin 
P(cosa)e  da  + p e , 


(2.33) 


This  equation  implies  that  the  pressure  is  continuous.  However,  for 
the  velocity  at  - 0 one  finds  from  Equations  (2.24)  and  (2  25)  that 
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P 6 

P(coscx)sina  -ikrcosi,  . *^o 
— ' e da  + 

pc 


ik.rcos<t>o 


P(.cosa)sina  -ikrcosa.  . *^o  o ik.rcos4)r.  .. .. 

— ^ -r e du  + — e , (2.J4) 

pc  pc 


which,  when  simplified,  becomes: 


n/  \ j -ikrcosu.  / , -Ikrcosu. 

P(cosa)sina  e da  = j -P(cosa)sinu  e du.(2.35) 


The  three  integral  Equations  (2-31),  (2.32)  and  (2.35)  are  a 
set  of  Integral  equations  that  satisfy  all  the  requirements  for  a 
solution  to  this  diffraction  problem.  Specifically,  these  equations 
satisfy  the  wave  equation  and  the  impedance  boundary  conditions  for 
an  incident  plane  wave.  To  solve  the  problem  one  must  determine 
some  function  P(cosa),  still  unkonwn,  that  will  satisfy  all  of  these 
integral  equations. 


2.5  The  Solution  of  the  Integral  Equations  for  Plane  Wave  Incidence 

Solving  for  the  function  P(cosa)  is  undoubtedly  the  most  onerous 
detail  of  this  study.  One  attempts  to  construct  the  unknown  function 
P(cosu)  by  considering  functions  with  properties  that  are  sufficient 
to  satisfy  each  of  the  integral  equations. 

Toward  this  end,  one  seeks  to  close  some  Y~path  within  a region 
of  convergence  and  then  evoke  the  Cauchy  Integral  Theorem  to  determine 
the  nature  of  P(co8a).  If  the  Cauchy  Integral  Theorem  is  to  apply, 
then  certain  restrictions  are  subsequently  Imposed  upon  the  functional 
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behavior  of  P(cosu).  From  these  restrictions  one  deduces  a functional 
form  for  P(cosa).  If  this  procedure  can  be  accomplished,  one  has  a 
solution  to  the  diffraction  problem. 

To  initiate  this  procedure,  one  must  specify  a y-path  and  the 
regions  of  convergence  through  which  the  y-path  may  be  closed. 

Starting  with  the  Integral  Equations  (2.31)  and  (2.32)  corresponding 
to  the  boundary  conditions  at  • ±v,  and  recalling  the  earlier 
section  defining  the  regions  of  convergence  in  Equations  (2.21)  and 
(2.22),  it  follows  that  these  regions  are 
for  < 0 

2nn  ^ (2  +1)tt  , 

and  for  > 0 

(2n-l)7r  ^ ^ 2nff  . 

Figure  2.5  shows  representative  y-paths  and  how  they  may  be  closed 
For  example,  one  such  path  originates  at  -loo,  travels  to  0,  then  to 
TT  and  terminates  at  n-ioo.  It  is  then  closed  at  infinity  along  the 
line  from  tt-Ioo  to  -ioo.  In  order  to  close  the  path  in  this  manner, 
it  will  be  convenient  if  the  Integral  along  the  part  of  the  path  at 
Infinity  is  zero.  For  this  condition  to  be  valid,  certain  restric- 
tions on  the  behavior  of  the  function  P(co8a)  at  infinity  must  be 
Imposed.  In  particular,  this  condition  maintains  that  the  absolute 
value  of  P(cosa)  must  be  bounded  at  infinity.  The  condition  on 


P(co8a)  is 
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P (cosa) 


M 


(cosa) 


as  ® 


(2.36) 


and  for  6 > 0. 

Now  consider,  for  example,  the  Integral  Equation  (2.31): 


L 


nt  \/  j \ Ikrcosa.  , j . 4,  X -Ikrcosu 

P(cosu)(sln0  +slna)e  dot  - p (sln0  -slm^  )e 

o o 


If  the  Y~path  Is  closed  In  the  manner  described,  one  finds  that  the 

Integral  around  a closed  path  equals  the  value  of  the  right  hand 

side  of  the  equation  which  Is  p^  (sln0^-sln0^)e  _ jjy  virtue 

of  the  Cauchy  Integral  Formula,  this  means  that  there  must  be  simple 

p 

poles  at  (2n-l)TT±(})Q  and  the  residues  are  -^^slnO^-slniJole 
for  some  y-path,  see  Figure  2.5.  Likewise,  the  Integral  In  Equation 
(2.32)  also  suggest  the  existence  of  poles,  but  with  a residue  of 
■^:j^(sln6  +sln(|)Q)e  ^ some  ypaths.  There  are  two  major 

features  about  the  behavior  of  P(cosa).  They  are: 

1.  P(cosa)  Is  bounded  for  . 

2.  P (cosa)  (sln0~±slna)  has  simple  poles  at  (2n-l)TT±(|)Q  and 

corresponding  residues  of  -:r-T(sin0”  sin4>  )e 

2TTi  o 

On  the  basis  of  this  limited  knowledge,  the  function  P(cosa) 
must  be  deduced. 

It  Is  difficult  to  proceed  without  giving  some  considerations 
to  the  physics  of  the  problem.  In  particular,  each  of  the  poles  at 
(20-1)711(1)^  have  distinct  physical  significance.  As  has  been 


^ to- 
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previously  discussed,  there  is  a multiplicity  behavior  in  the  complex 
a-plane,  so  the  poles  ^"*^0  ® family 

that  will  subsequently  be  repeated.  These  poles  represent  the 
following  physical  significance: 

a.  -TT-ip  represents  the  angle  of  the  specular  reflection 

from  the  lower  surface. 

b.  represents  the  angle  of  the  geometric  shadow  boundary. 

c.  Tr-0  represents  the  angle  of  the  specular  reflection  from 

the  upper  surface. 

d.  represents  the  angle  of  the  geometric  shadow  boundary 
For  example,  when  the  pole  of  Equation  (2.37a)  representing  the 
reflection  from  the  lower  surface  is  enclosed  by  a ypath.  Equation 
(2.32)  must  be  satisfied.  Likewise,  when  the  pole  in  Equation  (2  37c) 
representing  the  reflection  from  the  upper  surface  is  enclosed,  it 
must  satisfy  the  ether  Integral  Equation  (2.31).  To  accommodate  all 
these  requirements,  the  solution  is  considered  as  a product  of  two 
parts.  One  part  contains  the  necessary  singularities  while  the  other 
part  contains  an  analytic  function  that  adjusts  the  residues  to  the 
proper  values. 

Two  functions  that  have  the  proper  singularities  are  chosen  as 
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and 


sin 


g-TT 

2 


1 


sln4> 

o 

2 


(2.39) 


If  Che  first  function  in  Equation  (2.38)  is  enclosed  by  a Y~paCh 
that  includes  a simple  pole  at  the  reflection  from  the  upper 

surface,  one  must  modify  the  function  to  yield  Che  appropriate 
residue.  Such  a function  Chat  will  provide  the  proper  residue  is 
found  by  trial  and  error  to  have  the  form: 
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p cos<{i  (sin0  -sinS"*^)  (sin0  +sln<ti  ) ..  . 

*^o  o/2 _o o ^-lkrcos0Q 

2TTicos(t)^^2  (sinQ  +sin(})  ) (sin9'*'+sln(ti  ) 

o o 


p (sin4)  -sln9^)  . 

o ^-Ikrcosijio 

(sinij)  +sin0'*') 
o 


This  happens  to  be  exactly  what  is  necessary  to  satisfy  Equation 
(2.31). 

Consider  the  function  in  Equation  (2.39)  and  enclose  the  simple 
pole  at  i|)— •[r-(})^  for  the  reflection  from  the  lower  surface.  A suitable 
choice  that  will  satisfy  Equation  (2.32)  is 


p COS(]> 

o o 
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t 

I 


t 


I 


which  is  what  is  required  by  Equation  (2.32).  The  sum  of  the  two 

functions  in  Equations  (2,40)  and  (2.41)  then  represents  P(cosa) 

Separately,  the  two  functions  represent  the  contribution  “o  the 

diffracted  field  due  to  the  influence  of  the  upper  and  lower  surfaces, 

respectively.  In  the  parlance  of  Rlemann  sheets,  they  correspond  to 

the  solution  in  two  separate  sheets.  The  continuity  Equation  (2- 35) 

Implies  that  the  solution  P(cosa)  is  continuous  where  they  Join  or 

overlap;  namely,  at  the  shadow  boundaries.  As  a result,  if  the 

integral  is  closed  by  a Y”Path  that  encloses  the  geometric  shadow 

boundary  located  at  -ir+cj)  and,  again,  if  the  integral  is  closed  by  a 

o 

Y-path  surrounding  the  shadow  boundary  located  at  these  two 

roust  satisfy  Equation  (2,35). 

If  one  checks  this  requirement  for  the  assumed  func*'lons  of 
Equations  (2.40)  and  (2.41),  one  finds  that  the  residues  at  the  poles 
-TT+0^  and  in  Equation  (2.35;  are,  respectively: 


P cos$  ,,  . 

° ^ (sin0''^+8in4)  )(sin9  +sin<{i  )sin(;)  e o 

2 o o o 

2TTico80^  (slnO  +8ln(|)  ) (8ln6'''+8in(ti  ) 


p 8in<i 

o o -IkrcosipQ 

2ti1  ® 
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and 


p cosit  (sln0  -slnA  )(sln6  -sinit  )sin4< 
o o o o o 


^-lk.rcosi|)Q 


2i’icos0  (sin9  -sin({>  )(sin0  -sln^i  ) 
o o o 

2 

p sini}'  ^ ’ 

’^o  o -ikrcos<p_ 

• ■ e ^ , 

2nl 

From  this  result.  Equation  (2.33)  Is  also  satisfied.  Therefore, 
the  solution  for  P(cosa)  is  given  by 


p cos(i 
o o 
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All  that  remains  to  be  done  for  the  solution  of  the  diffracted  field 
Is  to  evaluate  the  Integral  In  Equation  (2.13)  using  the  solution 
shown  in  Equation  (2.42)  for  the  function  P(cosa). 

2 , 6 The  Far-Fleld  Solution  by  the  Method  of  Steepest  Descent 

To  find  the  far-fleld  solution  for  the  diffracted  pressure  when 
kr  is  large,  one  may  use  the  method  of  steepest  descent  which  is 
also  known  as  the  saddle  point  method.  Essentially,  the  method  is 
one  of  finding  the  saddle  point  and  Integrating  on  the  path  of 
steepest  descent  through  the  saddle  point.  Along  this  path,  the 
major  contribution  to  the  integral  will  result  near  the  saddle  point. 
If  the  y-path  is  closed  along  the  steepest  descent  path,  then  the 
value  of  the  integral  on  Y can  be  replaced  by  the  steepest  descent 
approximation  at  the  saddle  point  plus  any  enclosed  residues. 

Consider  the  integral  in  Equation  (2.13): 


p(r,ij>) 


V -ikrcos  ((f+ot) 
P(co8a)  e da 


36 


for  the  case  when  0<({)<iT  so  one  may  pick  the  upper  sign.  The 
point  Is  that  point  where  the  exponent  has  an  extremum,  l.e. 


’ (-ikrcos ((Ji-a) ) - 0 


-lkr8ln(0-a)  ■ 0 

and  the  saddle  points  are  located  at 

a “ i|'±2nTT  , 


The  corresponding  steepest  descent  paths  are  determined  by 


lm[-ikrcos  ) - 1 [-ikr]  = constant 


which  is 

-ikrcos (<()-u^)co8ha^  • -ikr 
or 

cos (0-a^)cosha^  - 1 


Write  Equation  (2.43)  as 


co8(4)-a  ) “ r — 

r cosha. 


saddle 

when 


(2.44) 


(2.45) 
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then,  for  -♦  , 

cos(({)-a^)  - 0 
or 

" i}>±(2n-l)  . (2  46) 

Equation  (2.46)  predicts  the  asymptotes  of  the  paths  as  they  tend 
to  Infinity. 

Finding  the  slope  of  the  path  Is  also  useful  In  locating  this 
path.  The  slope  Is  determined  from  the  following: 

2 

da^  sln((t>-a^)  cosh 

da  slnha^ 

r 1 


“ ± cosha^  . (2.47) 


From  Equation  (2.47),  the  slope  at  the  saddle  point,  ~ 0,  Is  1 
and  the  slope  of  the  asymptotes  are  As  a result  of  Equations 
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(2.46)  and  (2.47),  the  steepest  descent  path  is  shown  in  Figure  2.6. 
For  an  observer  angle  , the  steepest  descent  path  is  chosen 

according  to  Equations  (2.46)  and  (2.47)  and  the  requirement  that  it 
also  must  lie  in  a shaded  region  of  convergence.  It  should  be  noted 
that  the  pole  at  is  contained  within  the  closed  path. 

To  determine  the  far-field  solution,  use  the  Residue  Theorem  to 
evaluate  the  integral  over  the  Y~path. 

Using  the  Residue  Theorem  yields 


/ 


/ 

SDP 


, -lkrcos<J)-ci , ^ . -ik.rcos((t)-a)  , 

P(cosa)  e da  + / P(cosa)  e da 


- ni  Res(iT-$  ) « -2Tri  Res(TT-(l)  ) 
o o 


Then  using  the  steepest  descent  approximation,  Erdelyi  [14],  it 
follows  that 


P(cosa)  e *^^da  - - iri  Res(7r-4)  ) 

o 


-L 


. -ikrcos  (ifi-a)  , 
P(co8a)  e da 


p / sin<J)  -sin0^\ 


2 \ sim)!  +sin0 

o 


SDP 

lkrcos((|»-())o) 


+ 


kr  ! 47Ti  ® 


p cos 


0/2 


Figure  2.6.  Steepest  descent  path  for  the  far-field  pre 

7I-<()  < 4)  < TT. 
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for  kr  large 


. -ikrcos ((t>-a)  , *^0  „ -ikrcos ) 

P(cosa)  e da  - ~ e ^ 


I 


ITT,/ 4 -ikr  , 

e e p^cos4)o/2  ^ 

2>/  2TTkr 


sln4i+TT  - sln<}> 
2 


^-TT 

sin  - 


sin<i) 


(2  48) 


where 


sin6  - 8ln(l) 


sin4)  sin0 
o 


is  the  reflection  coefficient  and  'f.  are  defined  in  Equation  (2,43). 

1 , ^ 

It  is  interesting  to  note  that  the  residue  of  the  pole  that  is 
surrounded  by  the  closed  path  contributes  the  geometrical  optic  term, 
namely,  the  specular  reflection,  while  the  steepest  descent  path 
contributes  the  diffraction  term.  Consequently,  the  far-field 
diffraction,  without  the  geometrical  optic  terms,  can  be  expressed 
simply  as 


f 


171/4  -ikr 


Pj(r,({.) 


- « e PoCos4»o/2 


V2(4>) 


2/  2iTkr 


j8ln('^^]-  sini)) 


isin(-^^ 


■ sincf) 

o 

2 

(2.49) 
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Next,  consider  the  case  when  the  observer  angle  ^ is  0<((i<TT-(t  The 

o 

path  Is  closed  along  the  steepest  descent  path  but  the  pole  at  tt-j, 

o 

Is  surrounded  In  a different  manner  as  Indicated  in  Figure  2-7  This 
causes  a change  In  sign  of  the  geometric  term  la  Equation  (2.49)  and 
the  solution  In  this  region  becomes 


P(r  ,1)) 


[ n/  \ -ikrcos ((fi-a)  , 

I P(cosa)  e da 


C 

r 


^+ikrcos((li-(lio) 


+ P^(r.(^)  . 


(2  50) 


This  jump  in  the  specular  term  physically  represents  the  sudden 
transition  across  the  reflection  shadow  boundary  at  In  short, 

one  finds  that  there  are  three  distinct  regions  to  the  solution;  these 
regions  are  listed  in  Table  22  and  the  corresponding  physical  regions 
are  shown  In  Figure  2.8. 


Table  2.2.  The  three  regions  for  the  far-field  solution. 


Region  I 

IT  - 6 < 0 < TT 

O 

Region  11 

-TT  + A <(fc<ir-(b 
o o 

Region  111 
-V  < Ip  < -V  + (P 


^diffraction  ^ reflection  ^inc 
*^dif fraction  ^Inc 

^diffraction 
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The  far-fleld  solution  Is  helpful  in  searching  for  the  nulls  in 
the  diffracted  and  backscattered  field.  If  the  solution  is  considered 
in  the  form  shown  in  Equation  (2.49),  an  expression  for  the  zeroes  of 
the  diffracted  and  backscattered  pressures  may  be  found. 

For  the  case  of  diffraction  an  expression  that  locates  the  nulls 
may  be  determined  if  both  surfaces  have  the  same  real  impedance. 

This  expression  is  given  by 


sin())/  2 


_1  sln9 

2 sln(t^^2 


(2  51) 


where  0*0'*'-0  -pc/z. 

It  may  be  determined  from  Equation  (2.51),  that  when  the  impedance  is 
real  there  can  be  at  most  one  null  in  the  diffraction  function.  An 
expression  may  also  be  found  for  the  nulls  of  the  backscattered 
pressure.  In  this  case  there  can  be  at  most  two  nulls  in  the  solution, 
if  the  impedance  is  real.  When  the  impedances  are  the  same  on  both 
surfaces,  the  expression  for  the  location  of  the  nulls  is  given  by 


sinO 


sin0  . 


(2.52) 


2. 7 The  Near-Field  Solution  for  the  Diffraction  of  a Plane  Wave 
In  the  far-fleld  approximation,  one  notices  a Jump  at  the 
geometrical  shadow  boundaries.  In  reality  the  transition  across  the 
shadow  boundaries  is  continuous.  This  contrldlctlon  occurs  because 
the  steepest  descent  approximation  is  not  valid  if  the  saddle  point 
is  in  the  vicinity  of  a pole.  As  a result,  an  exact  solution  is 
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t 


needed  to  evaluate  the  solution  there  and  in  the  near-field  of  the 
edge 

To  achieve  an  exact  solution,  consider  P(cosa)  in  the  form 


P (coscj) 


cos 


(2.53) 


where  the  functions  <li  (a)  are  defined  as: 


4,  (a)  - / Aiiliill  + AiSlii).  + 

- T AttI  \ + 7 

^ ^ B(0’*a)  B(e^,a) 

, 


r.  « ^ 
A, 2 


^ A(9~.a) 
B(9'*',a)  B(9~,a) 


cos 


a±9 


^ C(9  .9'.a)  ' 

B(0‘^,a)*B(9“,a)  j 


(2.54) 


where 


A(9  ,a)  • (sin9  + slna) 


sln0  + 2cosa/2 


cos 


a-H)  ±tt1  :a-i)  nT\ 

o o 

— : ± sin, — : 


(2. 55) 


6(9”, a)  • (sin9“+  2co8a/2  cos<()  /2) 

o 


± sin 


sin9^+  2cosa/2[cos 


jO.+<t>  ITT 


±7T 

0 

. 

(2.56) 


and 


C(9^,9  ,a)  “ (8in9  -8in9  ) 8in(|)^*  (D*8ln9^*8in9  +[sln9^+sin9  ] 


•E  + 2co8a/2  co8(()  ;~*E)  , 
0/2 
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with 


cos 


a+<t)  +v 
o 


+ cos 


a+4)  -Ti 
o 


+ sin 


a-(t>  -Hr 
o 


a-4)  -TT 

- 4^|- 


and 


E * 2* cosa/ 2* 


cos 


I a-Hj)  +TT 

o 


+ sin 


a-(})  -htI  1 
o 


cos 


a+<t)  -TT 
o 


4 / 


- sin 


a-(t  -TT 

o 


4 

(2.57) 


A special  case  exists  when  one  surface  is  rigid;  i e. , 9 or  9 
equals  zero.  In  this  event  then,  and  can  be  simplified  into 
the  form 


= 4^ 


± (sin9'*'  sina)  (sin9'*'-t-2cosa/2*cos<J)^y2*®i''| — ^ ° 


a±TT±()) 


|sin6  +4*cosa/2*cosot/4"cos 

, d) 
atTT  + ^o 

cos 

1.0 


-TT 


sin9  -t-4cosa/2  cosq/4*cos 


sln6  +4cosu/2*cosa/4*cos 


4>/7T\1 


isin9  +2cosa/2*cos4) 

L 


(2.58) 


when  0 *0 


I 


I 
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and 


^2,1 (^) 


P 

o 

4TTi 


(sin6  Tslna)  sin9  +2cosa/2’cos0  ..•cos 

0/2 


a±TT+<ti 


♦ 


sin6  +4*cosa/2*cosa/4*co 


o/ 4 


-tt/4 


8in9  +4*cosa/2*cosa/4*cos 


1.0 


_^+TT^  a±7T+(})^ 

sin 


■1- 


sin9  +4*cosa/2*co8a/4*cos 

$ +TT\" 

^ 1 

sin9  +2*cosa/2*cos(}) 


0/2 
(2.59) 

when  9^  - 0°. 

The  singularities  of  P(cos  ) are  contained  in  the  terms  that 

appear  in  the  denominator  of  Equation  (2.53)  while  the  functions 

. .(u)  are  smooth  and  well  behaved. 

A t ^ 

One  seeks  to  find  the  near-field  solution  by  closing  the  Y~path 
along  the  steepest  descent  path  and  then  evaluating  the  integral  along 
the  entire  length  of  this  path  rather  than  Just  at  the  saddle  point. 

To  achieve  this,  consider  the  integral  along  the  steepest  descent  path 
as 


4 


48 


P (cosa) 


-Ikrcos (0-0)  , 
e do 


cos 


a+0ol 

2 ] 

1 COS 

/a-0. 


-ikrcos (0-a)  , 
e da  . 


(2.60) 


It  may  be  argued  that  since  0j^(a)  and  <I'2(a)  are  smooth  and  slowly 
varying  functions  near  the  saddle  point,  then  they  may  be  removed 
from  under  the  integral  in  Equation  (2.56)  and  evaluated  at  the 
saddle  point  a~0.  This  results  in  the  following  integral  along  the 
steepest  descent  path; 


■^DP 


P(cosa) 


-ikrcos (0-a) 


-ikrcos (0-a) . 
e da 


'SDP  cos  — 


'SDP  cos 


a-0 


(2.61) 


For  the  moment,  consider  only  the  first  integral  with  a change  of 
variable  a“0-a.  This  change  of  variable  will  affect  the  Y“Path,  the 
regions  of  convergence,  and  the  steepest  descent  path  merely  by  a 
lateral  shift  of  the  angle  0.  In  particular,  the  equation  for  the 
steepest  descent  path  becomes 


cosa 

r 


1 

cos  ha 

1 


(2.62) 


49 


This  change  of  variable  transforms  the  first  term  in  Equation  (2.61) 
to 


-ikrcosa 

e 


a-((>  -<i\ 

COS 

2 ) 

r 1 

1 

, 1 1 

COS 

a-'t'-'^oi 

o 

'sDP 

2 1 

cos , 

-ikrcosa  , 
e da 


(2.63) 


Using  the  steepest  descent  path  from  Equation  (2.62)  and  rearranging 
the  integrand.  Equation  (2.63)  becomes 


-ikrcosa 


da 


" <l>j^((}))*2*cos 


SDP 


-ikr-krsinha< sinur 
cosa/2  e ^ 

[cosa  + cos  ((})+(})  ) ] 
o 


da 


(2.64) 


It  also  follows  from  Equation  (2.62)  that 


sina  ■ tanhUj  . (2.65) 

r i 


Substituting  from  Equation  (2.65)  into  the  integral  in  Equation  (2.64) 


it  follows  that 
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ikrcosa  , 
e da 


cos 


a-<})-(t) 

o 

2 


4’j^(^)*2'cos 


-ikr  -krsinhaitanhai 
cosa/2'e  e ■ 


SDP 


[cosa  + cos(i})+4i^)  ] 


da 


(2.66) 


Since  the  exponent  in  the  left  side  of  Equation  (2.66)  decreases 
monotonlcally  away  from  its  maximum  value  of  1 at  a^  - 0 to  it 
is  convenient  to  change  the  variable  from  a to  T, 
where 

2 

slnha  tanha.  ■ T (2.67) 

o i 

such  that,  as  T runs  from  to  it  traverses  the  steepest  descent 
path.  Thus  the  steepest  descent  Integral  may  then  be  made  into  an 
Integral  whose  limits  extend  from  to  +<”. 

The  change  of  variable  that  is  needed  is  then 


t»<^2  sina/2 


(2  68) 


with 

. /"I  -iTT/4  , 

dx  “ — y-  e cosa/ 2 da  . 


(2.69) 


Substituting  Equations  (2.68)  and  (2.69)  into  Equation  (2  67)  it 


follows  that 
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-ikrcosa. 

I e da 

a-4.-4.  1 

SOP  cos  ^ — I 


<)>+<}>  ' /T  At 

. ...  . ^^0  1 cosa/2  e e 

41  (4i)  2 cos  — ^ I ~ 

^ ^ / "•  iTT  / 4 2 

[cosa+  cos  (4)+4'q)  ]co8<^/2  e -y 


-ikr  iTT/4^  , , X 2 _ . _ ol 

-e  e 4>  (4i)  cos  — — | 

1 ^ I 


dT 

9 4>+4>„  2 

[2cos  — — - 2sin  a/2]  , 


2 

,,  4)+4)  f -krt  , 

4,  cos  ^ ^ 

1 ^ ; 2 , ' 

j T +i‘2*co8  1— 2 — 


4>  (4>)e“^‘^’^e“^^^^2-b  / , ^ '2"  « 

^ j-co  [T^+ib^] 


(2.70) 


where 


b ■ /~2 


(2.71) 


The  Integral  in  Equation  (2.70)  can  be  expressed  in  terms  of 
the  complex  Fresnel  Integral.  Consider  the  integral 


1 N V-2 
J.»  f 


(2.72) 


I 
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with  V - kr.  From  Equation  (2-72),  it  follows  that 


d L -Ib^V 
dv  \ I 


■ -L 


-be 


-ib^v 


.,2 


dx 


/~ 


with 


t’ 


T 


(2.73) 


From  Equation  (2.73)  the  integral  in  Equation  (2.72)  may  be  written 
as 


I 


/"tT  be 


lb  V/  e 
V 


2 

-lb  V 


dv  . 


(2.74) 


2 2 

Letting  v ■ x /b  and  considering  b positive,  then  Equation  (2,74) 
becomes 


I - 2/T 


2Tr  F 


1 


ib  V -ix 


e dx  , 


U-HD 


cos 


-)  as). 


(2.75) 


where  F(a)  is  the  complex  Fresnel  Integral  defined  by 
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.2/  ,2 

F(a)  - e ^ / e ^ dx  . 


Some  useful  properties  of  F(a)  are  quoted  from  Clemmow  [A] 


iT  2 

F(a)  - e^^  / e dx 

s 


and 


r (a) 

F f I H 

F^(a)  - e / ® • 

J o 


the  following  identities  are  found  to  be  useful: 


F(a)  + F (a)  - i /T 

O 2 


and 


F(a)  + F(-a)  « /"iT  e ^ e^^ 


Since  the  Integral  in  Equation  (2.72)  is  an  odd  function  of 
corresponding  result  when  b is  negative  becomes 


I - -2  /T  F(-b/Tr  ) , 
so  consequently. 


(2.76) 

For 

(2.77) 

(2.78) 

(2.79) 

b , the 

(2.80) 


(2.81) 
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Substituting  the  expression  in  Equation  (2.81)  into  Equation  (2  66) 
and  finding  a similar  result  for  the  second  term  of  the  integral  in 
Equation  (2  bl),  it  follows  that  an  exact  expression  for  the  integra- 
tion along  C.he  steepest  descent  path  is 


P(cosa)e  . +4$  ((())/  it  e ^^’'e  2kr  cos4)+i}>  1 

1 I 2. 

I 2 

+4<l'.,(4i)/ TT  e~^‘^’^e"^^'^^Fl±/“2rf  cos({>-4)  i . (2  82) 

2 I o I 


i 2 

Again,  consider  the  first  term  of  this  expression.  Using  the  identity 
found  in  Equation  (2.78),  it  follows  that 


- +4*<1>^(()))/'‘ti 


±/~2kr  cos 

( 


l(P+(p 


/ TT  -iTT/4  i2k.rcos'‘ 
2 e e 


12krcos 

-e 


2 


1 

[0+4 

[ ±/  2kr  cos 

° 

2 

/ A iA.\ r~'  -ikr  -i7r/4  / tt  -iTT/4  ikr  ikrcos  (4i+(j)  ) 

+4*'1>^(<J))/  TT  e e ^ e e ^o 

_^ikr^ikrcos(<(i+(})Q)  C(±/  2kr  cos(4rt-())Q) 

-iS(i/~ 2kr  cos((|i+4)  )) 
o 


I 

1 ■ 


(2.83) 
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The  near-fleld  solution  for  Che  diffraction  may  now  be  written  as 
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/ /j.\  ikrcos ((ti+<ti  ) 

p^(r,4))  - ±2iTi$^((J))e  o 


1.0  - /T 


Itt/A 


<))+0 


C ±/  2kr  cos  ^ ■°  I -iS  ±/  2kr  cos 


0+41 


^1.0  - ^ 


±/1kr  COS 


<p+(p 


— -is  ±/~2kr  cos 


’ 1' 

1 2 

1 I . 

/ 

(2.88) 


2 . 8 Diffraction  of  Line  Source  Radiation 

Up  to  this  point,  the  type  of  excitation  has  been  a plane  wave. 
In  this  section,  a line  source  disturbance  will  be  considered.  The 
incident  radiation  from  a line  source  is  given  by 


Pi 


-If-/ A 
e 


P H 
o o 


(2) 


(kR)  -X-  p^ 


(2.89) 


for  kR>>l  and  where  R is  the  distance  from  the  location  of  the  source 

o o 

It  is  helpful  to  express  the  line  source  in  the  same  Integral 

form  of  the  preceding  sections.  Starting  with  a line  source  located 

at  X , y and  into  the  z direction  the  line  source  must  satisfy  the 
o o 


following  wave  equation: 
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+ k 


Pi 


4TTp  6(x-x  )6(y-y  ) . 
o o 


(2  90) 


Applying  Fourier  Transforms  on  x and  y,  it  follows  that 


- -0 


4it 


,222 

k -Y  -n 


where 


^00  oo 


* tt) 


^-in(x-Xo)^-iY(y-yo) 

“2  2 2 
k -Y  -n 


dYdn 


(2.91) 


Let  V “ 7^  and  B - kt/  1-v  , 

k 


then 


1 -ikv(x-Xo) 


Pi  “ If  Po  ; " 


I g-iY(y-yo) 

I 2 T~ 
J-oo  (e  -y") 


dY 


dv 


(2.92) 


The  inner  Integral  in  Equation  (2.92)  may  be  evaluated  using 
the  residue  theorem,  closing  the  path  of  Integration  from  -«  to  » 
in  the  lower  half-plane  and  enclosing  the  singularity  at  +6.  This 
yields 


P 


i 


ivk (x-Xq ) ^- iB (y-y o ) 
16 


kdv  . 


(2.93) 
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Next,  let  x“rcos(|),  x^“r^cos4)^,  y-rsinij),  y^-r^co3()>^ , and  v-cosa 
which  changes  the  integral  in  Equation  (2.93)  into  polar  coordinates 
and  gives  the  following  result: 


p g+ikroCos((()o-a)^- 

*^0 


Ikrcos  (<()-a) 


da 


(2.94) 


It  can  be  shown  that  the  integral  in  Equation  (2.93)  is  equal  to 
(2) 

TTiH^  ■^(kr).  With  this  knowledge  and  the  help  of  Equation  (2.89), 
the  line  source  has  an  Integral  representation 


Po" 


-i7T/4 


2tt 


^+ikroCos((t)o-a)^- ikrcos  (<))  a) 


(2.95) 


It  is  worth  noting  that  the  representation  for  a line  source  is  a 
convolution  and  is  obtained  by  multiplying  the  plane  wave  by 
e and  integrating  with  respect  to  a. 

If  such  a procedure  is  also  applied  to  the  plane  wave  solution 
for  the  diffraction,  a line  source  solution  will  result.  This  can 
be  shown  by  considering  the  Integral  equations  analogous  to 
Equations  (2.31)  and  (2.32)  for  the  boundary  conditions  for  a line 
source.  They  become,  respectively. 
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f ^ ~ ^ Tr  / A I 

I I _ , ^ ikrcosa  i e -ikr^cos ((ti„-£i)  ,„i  , 

/ I P(cosa,cosS)  (sinG  +sina)e  ! e o o dB'dot 


y J y 


27r 


-iir/A 


(2.96) 

/ 27r  - ^ 


and 


I i P (cosa.cosB) (sinG  -sina) 
J Y -'y 


-ikrcosa 


-Itt/A 


^-ikrgcos  ((pg-B) 


dB 


2it 


da 


-iTT/A 

-o! / e'^*^’^o‘^°®^'^’o"®^(sinG~-sinB)e^*^’^‘^°®^dB  . (2.97) 


/“ 2Tr 


Y 


Closing  the  path  in  the  region  of  convergence  and  using  the  residue 
theorem,  it  is  evident  that  the  same  functions  found  in  Equation 
(2  A2)  satisfy  these  integral  equations.  As  a result,  the  solution 
IOC  the  line  source  can  be  written  as 


P(r,^) 


-iiT/A 


27r 


4>^(a) 


Y ^YL 
■ikrcos(<))-a) 


i g+B  g-B 

cosl— 2“  cos  — ^ 


-ikrQCOs(<|)o-B) 


dpda 


(2.98) 


where  41  (a)  and  <l>2(a)  are  given  by  Equations  (2.53)  through  (2.57). 

Evaluation  of  the  integral  in  Equation  (2.98)  is  much  the  same 
as  the  procedure  for  the  exact  solution  derived  in  Section  2.7. 
Since  the  functions  4>j^«i)  and  4'2(“^  smooth  functions,  they  are 
evaluated  at  the  saddle  point,  a-i|>>  and  removed  from  under  the 
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Integral  sign.  The  resulting  Integral  Is 
-i7T/4 


p(r,ij>)  - 


2tt 


^-IkrQCos  ((J)o-B)g-ikrcos  (4i-a) 


Y J Y 


cos 


g+B 

2 


dudB 


-iir/A  I I -Ikr^cos (())n-B)  -ikrcos(4)-a) 
^ Jyj  Y 


Using  the  change  of  variable  4>^-6“B'  and  , the  first  term 

of  the  expression  in  Equation  (2.99)  now  becomes 

-iTT/4 


^1  *• 


$^(0)6 


/Tn' 


-ikr„cosB  -ikrcosa 


Y^  Y 


COS 

0+0^-a-6| 

2 ) 

dadB  . (2  100) 


Multiplying  the  numerator  and  denominator  by  cos 


'0+0^-g+B 


, then, 


1-0 
1 1 


(0)e 


•i  /4 


27T 


0+0  -g+Bj 

-IRr^cosB  -ikrcosa 
e o e 

cos(  2 j 

/<)+<()  -a-6) 
1 0 

COSj 

' 0+0^-g+B 1 

cos  1 ^ 

l 2 J 

dBda  , 


0^(0)e 


-iTT/4 


2/17 


y'y 


cosj 

0+0^-g+B 

-ikrocoag  - 
e ® e 

-Ikrcosg  .Q  , 
dBdg 

2 

cos 

/ 0+0^ -g 

-8in6/2j 

cos 

fdH)  -a 
0 

+sinB/2] 

J 

1 2 

1 2 , 

2/17 


-iTr/4 


-a\  1 '*’  ■ -ot 

Y 9 sinB/2|  cos  ^ — +sin6/2 


-IkroCosB  -Ikrcosa 


dad  3 


(2.101) 


In  the  first  term  of  the  integrand  in  Equation  (2.101),  change  a to 
-a,  then 


4 r -y-^  cosa/2cos6/2  V^‘^""°"“d6da 


1 )j  2 11  r 

Yy  jcosa+cosS-2-4sin— ^ina/2sinS/2+2cos  ° 

(2.102) 


Now  use  the  change  of  variable 


C “ /I  e ^^^^8lna/2  and  n - / 2 e ^^^^8inB/2 


(2.103) 


which  results  in 


4$^(^)e-^^/''  -ik(r+r.) 

I.  - cos—: — e 

^ / 27  2 


2 2 

2 

>H>o 

2 

nf^'^icos 

2 

(2.104) 
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Using  another  change  of  variable  to  p and  y,  where 


/ h / h 

5 " V — pcosY  and  h ■ V “ 
o 


psiny 


(2.105) 


with 


R)  - (r+r^). 


It  then  follows  that, 


(2.106) 


^1- 


. -lTr/4 

— 4>.  ^1  cosl 

/~w  ^ ' 


l4H-<p 


pK(p)e'‘^lP  dp  , (2.107) 


where 


K(p)  - 


dY 


/ 2 
' P 


fx — 2 / 2 

7—  cos  y+J  — sin  Y+2sin 
r r 

o 


slnycosY 


[FT  M**): 


1/ 


(2.108) 


The  quantity  K(p)  was  evaluated  according  to  Clemmow  [4],  such  that 


Zn  sec 


K(P)  - 


(P-Hp 


/ o o 

p +21p cos 


(j)+4, 


(2.109) 
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Substituting  Equation  (2.109)  into  Equation  (2.107)  ylalds 


±4>^  (0) 4 • e" ^e“ 


J 


r“  ..-‘^ip 


£±. 


0^+1 

[Vsll 

p 

«JI 

,(2.110) 


/ 2 2 I “^"^o 

with  S - / r^+r  -2rr  cos  j--!  and  ± for  cos 
o o 1 2 


0 . 


2 2 

Set  iX  - -ik.(Rj^-S)  + kRj^p  , then  the  integral  can  be  put  into  the 
form 


I - It  ((>)*4‘<^2Tr  e ^ 

^ ^ ^ 2 

J ✓ X +2kS 


dX 


(2.111) 


When  k(r+r^)>l,  an  approximation  for  Equation  (2.111)  can  be  made 
by  replacing  the  non-exponential  part  of  the  integrand  by  the  value 
of  X at  the  lower  limit  and  removing  it  from  under  the  integral  sign, 
such  that 


^1-  ^ 


/ A fA.\  /~T“  -I't/A  -ikRi 
4*'l>,  (i^)/  2tt  e e 


e dX  , 


/ k(R^+S) 


4(j> 


g-iiT/4^-lkRi 


/ k(Rj^+S) 


± 


F(/  k(Rj^-S)  , 
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47Ti4>  / 1 

^ ;1.0  - / 2 e ' \ C(/  k(R  -S)  ) -ISC/TTr.-S)  ) , 

/■2krRj^+S)  1 I 1 

(2.112) 


A similar  expression  can  be  determined  for  the  second  term  1^ 
in  Equation  (2.92).  As  a result,  the  solution  for  the  line  source 
when  k(r+r  )>1  can  now  be  written: 


lyjy 


4^1$^  (4)) 

/ 2k(Rj^+S) 


1.0  - /~2 


■1  /4 


C(/' k(R^^  )-iS(A(R^-S) 


4TTi<I>^(4»)e 


-kR 


/ 2k(Rj^+R)  [ 


1.0  - /T 


i7T/4 


|c(/  k(R^-R)  -iS(/  k(R^-R)  )| 


+ (Geometrical  Optic  Terms) 


where 


R - /^^+r  ^-2rr 


cos((l)-<|)  ) 
o o o 


2^  2 , 
/ r +r  -d 


S - / r"+r  “-2rr  co8(4)+4)  ) 
0 0 o 


and 


R,  - r+r 
1 o 


(2.113) 
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An  asymptotic  approximation  for  the  line  source  solution  may  be 
made  if  the  location  of  the  line  source  and  the  observer  are  both  very 
far  from  the  half-plane,  i.e.,  kr  ->•  " and  kr^  Applying  the 

method  of  steepest  descent  to  the  double  integral  in  Equation  (2.89), 
the  far-field  approximation  becomes 


P (rA)  “ 


0 

(4>) 

$2(4>) 

2 )/~2'nkr  ^2t\V.t 

0 

cos 

f^-N-ol 

1 2 1 

1 

cos  . 

\ ^ 

(2.11A) 


which  may  also  be  written  as 


/ 2TTkr  /Tffkr 

o 


\ 

*'^(<J>) 

^2(*>  1 

sin 

1 

1-  sin  J 
2^ 

T 

;-r-ri 

!-  J 

(2.115) 


2.9  The  Diffraction  of  Point  Source  Radiation 

The  diffraction  caused  by  a point  source  excitation  can  be 
determined  by  considering  the  integral. 


p(r,()),z)  - 


P(cosa)e 


-ikt/  r^+ro^+2rrQCOs(<()-a)+(z-ZQ)^ 


— • (2.116) 


k/  r^+r  ^+2rr  cos(4)-a)+(z-z 
I o o o 


where  P(cosa)  is  defined  as  in  Equation  (2.A2).  By  analogy  to  the 
plane  wave  solution,  this  can  be  regarded  as  a spectrum  of  point 
sources.  A point  source  located  at  r^,  and  z^  is  expressed  as 
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.z) 


9 9 ' 

-ik/  r +r  “2rr  cos((t)-4)  )+(z-z  )' 
p e o o o o 

*^o 

k/  r^+r  ^-2rr  cos (0-0  )+(z-z 
0 0 o o 


(2.117) 


Applying  the  boundary  condition  on  both  surfaces  yields  integral 
equations  analogous  to  Equations  (2.31)  and  (2.32)  which  are  expressed 
respectively  as: 


j + r 

P(co8a)  ( sinO  +r  sina  1-  rrr  ; -r-r — da 
1 o kR  / kR 


r 

-p  {sin9  -r  sin0  1-  — ! ) ^ . 

0 1 o ^o  kR  : / kR* 


isinO  -r  sin0  [l-  tt"]  1 


r sin0  1-  — / — da 
0 0 kR  / kR 


P(cosa)  / 


^ ' -ikR* 

-p  (sinS  +r  sin0  1-  r-rr  ) — > 
o 0^0  kR*  / kR*  ’ 


(2.118) 


(2.119) 


where 


R " r^+r  ^-2rr  co8a+(z-z  )^ 


+r  -2rr  co8a+(z-z  ) , 

o o o ’ 


/ 2 2 2 
R*  - / r +r  -2rr  co80  +(z-z  ) 
o o o 


If  the  'y~path  is  closed  in  the  region  of  convergence  and  the 
residue  theorem  is  applied  to  the  singularities  enclosed,  then  the 
function  P(cosa)  of  Equation  (2.42)  will  indeed  satisfy  Equations 
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(2.118)  and  (2.119).  This  implies  that  the  result  in  Equation  (2.116) 
is  the  diffraction  solution  for  a point  source  excitation. 

Because  the  exponent  in  the  integrand  for  the  point  source  differs 
from  plane  wave  integral,  it  is  helpful  to  check  the  regions  of 
convergence.  Also,  because  of  the  square  root  term  in  the  denominator, 
branch  points  must  be  accounted  for. 

Consider  the  function  in  the  exponent 


R(a)  - / r^+r  ^+2rr  cos  (((i-a)  + (z-z  )^ 

O 0 0 


(2.120) 


The  branch  points  are  located  at  values  of  a where  R(a)  » 0,  which 
implies 


2 2 2 
r +r  +(z-z  ) 

0 0 


- cos(4i+TT-a) 


(2.121) 


where 


“ cosha^'cosa^'  , 


a'  - TT+(|)-a  . 


2 2 2 2 2 

Since  (r-r  ) - r +r  -2rr  > 0,  then  r +r  > 2rr  . It  follows  that 

O O 0—0 
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'22  2 
I r +r  +(z-z  ) 

I o o 

I 2rr 

I o 


> 1. 


(2.122) 


As  a result  of  Equations  (2.122)  and  (2.121),  the  branch  points  are 
In  the  complex  plane  and  located  at 

a ' - 2niT 
r 

and 

a^'  - ± cosh 

' o 


2^.  ,2^ 

r +r  +(.z-z) 
o o 


2rr 


±b 


or 


- (2n-l)Tr+<)) 


and 


- ±b 


(2.123) 


The  branch  cuts  originate  at  a “ (2n-l)TT-H|)±ib  and  extend  along 
straight  lines  to  infinity  at  a ■ (2n-l)iT+<()±i».  The  function  R(a) 
can  be  written  in  the  form 


R(a)  - |r|^^^  e^"^  - |r|^^^  (co3iJH-lsini|))  , 


(2.12A) 
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where 


2rr  slnha, sln(0-a  ) 

tanti»  - - — ^ . (2.125) 

r +r  +2rr  cosha.cos((ti-a  ) 
o o i r 

To  guarantee  convergence,  the  condition  sini^)  < 0 is  required.  This 
is  necessary  so  that  the  exponent  will  reflect  a decaying  behavior 
and,  thus,  allow  convergence  of  the  integral  in  Equation  (2.116). 

The  Table  2.3  helps  to  locate  the  regions  of  convergence  by 
listing  the  values  of  a such  that  R(a)  is  real.  These  values  will 
determine  the  boundary  between  the  regions  of  convergence  and  the 
regions  of  divergence.  Consequently,  from  Equations  (2.124),  (2.125), 
and  Table  2.3,  the  branch  points  and  regions  of  convergence  can  now 
be  Illustrated  in  Figure  2.9. 


Table  2.3.  The  values  of  ar  and  that  describe  the  curve 
in  the  a-plane  as  R(a)  traverses  the  real  axis. 


Or  R(a) 


(2n-l)  TT  + $ 
(2n-l)  IT  + (() 
(2n-l)  TT  + 4> 
(2n+l)  TT  + 4) 
(2n+l)  TT  + 4» 


+ 00  + 

+ b 0 

0 ± /(r+r  )^  + (Z-Z  )^ 

o ^ 
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Evaulatlon  of  the  Integral  In  Equation  (2.116)  proceeds  in  a 
manner  similar  to  that  for  the  plane  wave  and  the  line  source.  Using 
the  function  P(cosa)  in  the  form  of  Equation  (2.53),  the  smooth  part 
of  the  function  is  evaluated  at  the  saddle  point  and  removed  from 
under  the  Integral  sign  so  that 


Consider  the  fir  at  term  and  cose  the  yp^th  along  the  line  i”  +4)  to 
-1“  +ip.  This  results  in  the  following  integral  for  the  first  term: 


/2  2 2 
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2cos 
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k/  [ (r+r^)^+(z-z^)^] 


./i 


1- 


4rr  sinh  a. 
o 1/2 

(r+r^)^(z-z^)2 


(2.127) 


Using  the  change  of  variable  x 
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— - e ^ sinhaj^^2»  where 
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“ /(r+r^  +(z-z^)  , and  applying  the  Identities  cos 
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-1  and  cosha^  ••  l-28inh  (a^/2),  the  integral  in 


Equation  (2.127)  becomes 
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Expand  the  exponent  term,  kR^^/  1-x  , in  a Taylor  series.  If  kRj^>  1, 

retain  only  the  first  two  terms  of  the  expansion,  evaluate  the  slowly 

2 1/2 

varying  function  (1-ix  ) in  the  denominator  at  the  saddle  point 


J 
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T = 0,  and  remove  It  from  under  the  Integral  sign.  As  a result,  for 
> 1,  the  Integral  becomes 


l<t^(4i) 

kR, 


,,  I 2/krr 
+1it/4  I o 

e — :: cosl 


<P+(P 
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which  Is  In  the  form 
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The  integral  in  this  form  is  identical  to  that  of  Equation  (2.68) 
which  now  yields 
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There  is  a similar  expression  for  the  second  term;  thus,  the 
exact  solution  becomes 


P„(r.4.,z) 


TTl'I'^((|))e 


ikrr 

■ikR^  Tp 

e 


kR, 


X 


1.0-  /2  e 


i /4 


-iS 


(2.132) 


75 
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The  asymptotic  approximation  for  the  point  source  may  also  be 
made  when  kR^^  kr  ■+  »,  kr.  In  this  event,  Equation  (2.116)  may 

be  simplified  to 


lk2rr 


p(r,((),z) 


-Ikr, 


2 (t-a 
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(2.133) 


Again,  the  steepest  descent  point  Is  a - ()).  The  far-fleld  approxi- 
mation may  then  be  written  as 


p(r,  ,z)  - 


-IkR 

2Tr  P(cos»)e^  ^^e 


kR, 


kr  / kr. 


(2.134) 


where  P(cos4>)  can  be  found  In  Equation  (2.42). 


CHAPTER  III 


NUMERICAL  RESULTS 


3. 1 Introduction 

Numerical  results  for  various  source  type  and  location,  and 
Impedance  cover  conditions  are  presented  in  this  chapter.  The  para- 
meters are  varied  to  explore  the  dependence  of  the  diffracted  pressure 
and  the  backscattered  pressure  on  the  impedance  cover.  Calculations 
are  also  made  for  a variety  of  source  and  observer  locations.  Special 
attention  is  given  to  the  far-fleld  approximation.  All  graphs  are 
polar  plots  of  the  angular  distribution  of  the  diffracted  or  back- 
scattered  pressure  measured  in  decibels  relative  to  the  amplitude  of 
the  Incident  pressure. 

3.2  Impedance  ConBiderations 

In  the  solution  for  the  diffracted  pressure  for  all  three  source 
cases,  the  parameter  that  determines  the  impedance  condition  is  the 
variable  6^  or  6 . From  Equation  (2.3U),  one  recalls  that  sin(6)~pc/z. 
The  various  impedances  that  are  considered  in  this  study  and  their 
corresponding  values  of  6 are  listed  in  Table  3.1. 


Table  3.1.  Impedance  Conditions  and  their  corresponding  Brewster 
Angles  6. 


Z 

pc/Z 

6 Brewster  Angle 

Pressure  Release 

0 

00 

ioo 

Matched 

pc 

1.000 

o 

o 

l.lbpc 

.866 

o 

o 

2. OOpc 

.500 

O 

O 

Rigid 

oo 

0 

0° 
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CompuCatiuns  were  performed  using  Che  Impedances  listed  In 
Table  3.1  for  specific  reasons.  The  cases  when  6 - 51  or  6 • 0*^  are 
for  Che  classical  cases  of  pressure  release  (soft)  and  rigid  boundary 
conditions  respectively.  The  case  of  6 > 90^^  is  particularly 
Interesting  since  this  case  represents  a totally  absorbing  surface 
The  other  cases  are  computed  to  complete  a range  of  examples.  For 
each  of  the  different  source  cases,  the  far-fleld  pressure  has  the 
same  angular  dependence  as  given  by  Equation  (2.4^).  This  angular 
dependence  Is  calculated  to  Indicate  only  tlie  pattern  of  the  far-fleld 
pressure.  It  is  possible  to  have  at  most  one  zero  occurring  In  the 
far-fleld  diffracted  pressure  and,  at  must,  t«ro  nulls  occurring  In 
the  far-fleld  backseat tered  pressure  The  Influence  of  the  Impedance 
on  the  location  of  nulls  Is  also  denwnstrated  by  these  computations 

3. 3 Source  nbserver  Considerations 

The  location  of  the  obsetver  la  varied  In  an  attempt  to  determine 
the  effect  of  the  radial  parameter  and  the  angle  of  Incidence  on  the 
behavior  of  the  diffracted  or  backscattered  pressure  The  dimension- 
less variables  kr,  ki  , and  k(z-z  ) are  considered  for  a wide  range 

o o 

of  values  to  demonstrate  the  transition  from  the  near-field  to  the 
far-fleld.  Comparisons  are  also  made  between  the  three  types  of 
source  conditions  when  the  radial  variables  are  considered. 

In  the  case  of  diffraction,  the  angle  of  Incidence  was  chosen 
to  be  ■ 120°  as  a representative  case.  For  backscatter,  the  angle 
of  Incidence  and  Che  observer  angle  are  equal.  In  other  words 
<))  ~ In  Equation  (2.42)  for  Che  backscatter  pressure. 
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3. k Diffraction  and  Backscatter Ing  of  a Plane  Wave  by  an  Impedance 
Covered  Kalf-Plane 

All  of  the  graphs  for  the  diffracted  and  backscattered  pressure 

are  In  polar  form  and  are  measured  In  decibels  relative  to  the 

amplitude  of  the  Incident  pressure.  In  the  exact  solution  for  the 

plane  wave,  the  factor  > kr  appears  in  the  argument  of  Che  Fresnel 

Integrals  as  shown  by  tquaClor  (2  86).  Consequently,  Che  magnitude 

of  the  pressure  will  diminish  at  a rate  proportional  to  l//kr. 

because  of  this.  It  Is  desirable  Co  normalize  Che  pressure  by  Che 

same  factor  in  order  to  facilitate  the  comparison  of  successive 

diffraction  plots  As  a result,  what  is  really  being  sltown  In  Che 

plots  is  (r,t)*>'kr.  In  effect,  one  Is  observing  only  Che  angular 

dependence  of  the  diffracted  or  backscatCered  pressure.  Throughout 

this  -hapCer  t fie  term  "illuminated"  is  borrowed  from  the  field  of 

optics  and  corresponds  to  the  expression  "insonificd"  in  acoustics. 

There  are  fourteen  figures  that  deal  with  the  diffraction  of  a 

plane  wave  by  a fialf-plane.  In  all  of  these  plots,  the  angle  of 

incidence  is  set  at  ^ • 120^ 

o 

The  Influence  of  the  impedance  of  the  illuminated  surface  Is 
considered  in  the  first  eight  figures.  These  figures  are  plotted 
for  Che  following  parameters: 

1.  Figures  3.1  and  3.2  have  an  Impedance  on  the  Illuminated 
surface,  6^  ■ 30°,  and  6 * 0 (rigid)  on  Che  unllluminated 
surface.  In  Figure  3.1,  kr  is  20  and  50  and  In  Figure  3.2, 
kr  Is  100  and  infinity. 
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^ 2.  Figures  3.3  and  3.4  have  G"*"  = 60°  and  0 ”0°  (rigid). 

The  values  for  kr  are  20  and  50  in  Figure  3.3;  and  100,  and 
infinity  in  Figure  3.4. 

3.  The  impedance  conditions  are  G"*”  = 90°  (totally  absorbing) 
and  0 - 0°  (rigid)  in  Figures  3.5  and  3.6.  In  Figure  3.5, 
kr  is  20  and  50  and  in  Figure  3.6,  kr  is  100  and  infinity. 

4.  Figures  3.7  and  3.8  represent  the  situation  of  pressure 

release  on  the  illuminated  surface  G"^  = 15  and  rigid  on  the 
opposite  surface  0 =0°.  The  values  of  kr  are  20  and  50  in 

Figure  3.7,  and  100  in  infinity  in  Figure  3.8. 

The  next  four  plots  demonstrate  how  the  impedance  on  the 

unilluminated  surface  influences  the  diffracted  pressure.  The  angle 

of  Incidence  is  d*  = 120°. 

o 

Figures  3.9  and  3.10  represent  the  diffraction  by  a half-plane 
with  pressure  release  Impedance  on  the  unlllumlnated  surface,  l.e., 

G"*^  “ 5i,  while  the  illuminated  surface  has  0 =0°  (rigid).  Again, 
kr  - 20,  50  in  Figure  3.9  and  kr  - 100  and  infinity  in  Figure  3.10. 

The  Impedance  condition  on  the  upper  surface  in  Figures  3.11 
and  3.12  is  G"*"  ■ 90°  (totally  absorbing)  while  0 -0°  remains  the 

same.  In  Figure  3.11,  kr  - 20  and  50,  while  in  Figure  3.12,  kr  » 100 
and  infinity. 

The  Influence  of  an  imaginary  component  of  Impedance  is  con- 
sidered in  Figures  3.13  and  3.14  wtiere  G^  - 30  + i.l  and  0 • 0°. 

The  angle  of  incidence  is  ((i^  “ 120°  and  kr  - 20  and  50  in  Figure 
3.13,  while  kr  ~ 100  and  infinity  in  Figure  3.14. 
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raction  of  a plane  wave  with  6 *30  i.l.  0 • 0 . an( 
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The  next  ten  plote  are  concerned  with  Che  beckecetCered  preeeure. 
For  this  situation,  the  source  and  observer  are  at  the  sane  position, 

l.e.  it  > in  Equations  (2.42)  and  (2.32).  The  dependence  of  the 
backscactered  pressure  on  the  impedance  cover  can  be  studied  by 
comparing  Che  following  plots: 

1.  The  impedance  conditions  for  Figures  3.13  and  3.16  are  6^ 

- 30*^  and  6 - 0.  The  variable  kr  is  20  and  30  in  Figure  3.13, 

and  kr  is  100  and  infinity  in  Figure  3.16. 

2.  In  Figures  3.17  and  3.18,  6^  - 60°  and  6 ■ 90°  while  kr  ~ 20 

and  30  in  Figure  3.17,  and  kr  > 1(X)  and  infinity  in  Figure  3.18. 

3.  Figures  3.19  and  3.20  have  a totally  absorbing  upper  surface 

with  6^  - 90°  while  6 - 0°.  In  Figure  3.19,  kr  is  20  and  30, 

while  kr  is  100  and  infinity  in  Figure  3.20. 

4.  The  Impedance  conditions  in  Figures  3.21  and  3.22  represent 

a pressure  release  upper  surface  6^  > 13  and  a rigid  lower 
surface  6 ••  0°.  Again,  kr  is  20  and  30  in  Figure  3.21,  and 

kr  Is  100  and  infinity  in  Figure  3.22. 

3.  In  Figures  3.23  and  3.24,  the  impedance  cover  is  the  same 
on  both  sides  with  6^  • 6 > 30°.  The  values  of  kr  in  Figure 

3.23  are  20  and  30,  while  in  Figure  3.24  they  are  100  and 
infinity. 

3.5  Diffraction  and  Backscatterlng  of  Line  Source  Radiation  by  an 
Impedance  Covered  Half-Plane 

Computations  for  Che  line  source  are  made  primarily  to  study  the 
effect  of  Che  source  type  and  location  on  Che  diffracted  and  back- 
scattered  pressure.  In  order  Co  compare  these  graphs  with  Che 
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diffraction  of  a plane  wave,  the  angle  of  Incidence  is  kept  at 

4>  ••  120*^.  For  the  line  source  the  magnitude  of  Che  pressure 

o 

diminishes  at  a rate  proportional  to  l//kr*l//kr^  as  Implied  by  the 
argument  of  Fresnel  Integrals  in  Equation  (2.113)  and  the  demonlnator 
of  Equation  (2.114).  Consequently,  Che  calculacions  for  Che  line 
source  are  normalized  by  the  factor  /kr*/kr^  in  order  to  permit  a 
better  comparison  of  results. 

The  next  four  figures  are  concerned  with  Che  diffraction  of 

line  source  radiation.  They  are  shown  in  an  attempt  Co  illustrate 

the  influence  of  kr^  on  the  behavior  of  the  diffracted  pressure. 

The  values  of  kr  are  chosen  to  be  kr  • 20,  50,  100  in  order  to 
o o 

provide  a range  of  values. 

In  Figures  3.25  and  3.26,  Che  impedance  on  Che  two  surfaces 

are  O"*"  - 60°  and  9 "0.  The  value  of  kr  is  fixed  at  50,  while 

kr  - 20  and  50  in  Figure  3.25  and  kr  ~ 100  in  Figure  3.26. 
o o 

The  impedances  in  Figures  3.27  and  3.28  are  9^  ■ 90°  (totally 

absorbing)  and  9 >0.  Again,  in  each  of  these  figures,  kr  • 50. 

In  Figure  3.27,  kr^  - 20  and  50,  while  in  Figure  3.28  kr^  - 100. 

The  next  three  figures  represent  Che  backscatCered  pressure 

of  a line  source.  For  this  situation,  Che  source  and  observer  are 

at  Che  same  location;  l.e.,  kr  ■ kr  and  4>  ~ Equation  (2.113). 

o o 

In  each  of  the  following  figures,  Che  Impedances  on  Che  two  surfaces 

are  equal.  In  Figure  3.29,  the  impedance  conditions  are  9^  « 9 - 30°, 

while  kr  - kr^  - 20,  and  kr  - kr^  “ 100.  The  Impedances  in  Figure 

3.30  are  9^-9  - 60°.  Again,  the  radial  parameters  are  kr  - kr^ 

- 20,  and  kr  - kr^  - 100.  In  Figure  3.31,  the  impedance  conditions 

are  8^  ■ 9 ■ 90°,  while  kr  ■ kr  ■ 20  and  kr  ■ kr  ■ 100. 

o o 
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3. 6 Diffraction  and  Backscatter Ing  of  Point  Source  by  an  Impedance 
Covered  Half-Plane 

The  last  type  of  source  to  be  considered  Is  the  point  source. 

For  this  case.  It  Is  necessary  to  specify  the  relative  distance  In 

the  z-dlrectlon,  k(z-z^),  between  the  source  and  the  observer,  as 

well  as  the  parameters  kr  and  Itr^.  The  dependence  of  the  diffracted 

pressure  and  backscattered  pressure  on  k(z-z^)  Is  explored  In  this 

last  section.  The  magnitude  of  the  pressure  for  a point  source 

decays  at  a rate  proportional  to  ’l/Air  as  demonstrated 

1 o 

In  Equations  (2.133)  and  (2.135).  Consequently,  the  calculations 
are  normalized  by  this  factor  In  order  to  facilitate  comparisons 
between  the  different  plots. 

The  diffracted  pressure  Is  considered  In  the  next  group  of 

figures.  The  parameter  k(z-z^)  takes  on  the  values  0 and  100.  The 

purpose  of  using  these  values  Is  to  Investigate  Its  Influence. 

There  are  two  combinations  for  the  variables  kr  and  kr  that  are  used. 

o 

these  combinations  are  kr  - kr  - 20,  and  kr  - kr  - 100.  Again 

o o 

this  choice  Is  made  to  provide  a range  of  cases  suitable  to  explore 

their  Influence  on  the  behavior  of  the  diffracted  pressure.  The 

graphs  are  displayed  In  such  a manner  that  one  may  examine  the  trends 

as  k(z-z  ) varies  and  as  kr  and  kr  vary.  These  calculations  are 
o o 

then  made  for  two  representative  impedance  conditions.  In  both 

cases,  the  Impedance  on  the  lower  surface  Is  fixed  at  6 ~ 0*^,  while 

the  Impedance  on  the  Illuminated  side  has  the  values  6^  - 60°  and 

90°.  The  location  of  the  source  Is  fixed  at  ((i  ■ 120°. 

o 

The  diffracted  pressure  of  point  source  radiation  Is  Illustrated 
In  the  next  six  plots. 
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Figure  3.32  shows  the  diffraction  by  a surface  with  Impedance 

conditions  0^  • 60°  and  0 - 0°  for  kr  - kr^  ■ 20.  In  this  figure 

k(z-z  ) “ 0 and  100. 
o 

Figure  3.33  has  the  same  surface  condition  with  0^  - 60°  and 

0 *0°  but  kr  “ kr^  • 100.  Again,  one  may  Judge  the  Influence  of 

k(z-z  ) since  It  has  the  values  ol  0 and  100. 
o 

The  Impedance  conditions  In  Figure  3.34  are  once  again  O'**  - 60° 
and  0 “0°  but  k(z-z^)  - 0.  One  may  consider  the  effect  of  kr  and 

kr  since  they  take  the  values  of  kr  - kr  ■ 20  and  kr  - kr  ■ 100. 

O 0 0 

In  Figure  3.35,  0"*"  - 90°  (totalling  absorbing)  0 -0°  and 

kr  - kr^  - 20.  The  two  plots  In  this  figure  are  for  k(z“z^)  ■ 0, 

and  k(z-z  ) ■ 100.0. 
o 

Figure  3 36  shows  a half-plane  with  0"*"  » 90°  and  0 - 0°  with 

kr  - kr  - 100. 
o 

In  the  last  figure  for  diffraction,  O"*"  - 90°  and  0-0°  while 

k(z-z  ) “ 0.  The  plots  In  Figure  3.37  are  for  kr  ••  kr  - 20  and 
o o 

kr  - kr  ■ 100.  They  Illustrate  the  effect  of  kr  and  kr  . 
o o 

The  final  series  of  figures  Investigates  the  backscattered 

pressure  of  an  Impedance  covered  half-plane  for  a point  source 

excitation.  For  backscatter,  the  source  and  receiver  are  at  the 

same  location,  so  k(z-z^)  Is  always  zero.  Computations  are  made 

for  kr  - kr  > 20,  and  kr  ■ kr  " 100.  These  values  for  kr  and  kr 
o o o 

represent  the  cases  when  the  source  Is  close  to  the  edge  and  also 
when  It  Is  far  from  the  edge.  The  two  representative  Impedance 
conditions  used  for  this  group  of  figures  are  0^  0 - 30°,  and 

0"^  - 0"  - 60°. 


rR*=:  IDO  K<Z-Z*D^D  KR=  rRt=1DD  ICCZ-Z«D=1DD 


Figure  3.35.  Diffraction  of  point  source  radiation  with  0 
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Figure  3.38  shows  Che  backscattered  pressure  for  kr  ~ kr^  ~ 20 

and  kr  - kr  - 100  when  0^-0  - 30°. 

o 

Figure  3.39  whows  the  backscattered  pressure  for  kr  ■ kr^  • 20 
and  kr  • kr  - 100  when  0^-0  ••  60. 


CHAPTER  IV 


DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 


4 1 Introduction 

The  major  task  of  this  Investigation  Is  to  model  mathematically 
a certain  physical  situation  and  derive  a solution  based  upon  that 
sodel.  Specifically,  the  problem  of  the  diffraction  of  acoustic 
waves  by  a locally  reacting  impedance  covered  half-plane  is  solved. 
The  solution  resulting  from  this  analysis  Is  significant  in  that  It 
Is  a closed  form  solution.  Being  a closed  form  solution  facilitates 
the  actual  physical  interpretation  of  the  problem.  A discussion  of 
the  physics  of  the  diffraction  problem  solved  In  Chapter  II  follows 
In  Section  4.2.  Additionally,  some  specific  cases  are  examined  for 
the  Influence  of  the  Impedance  cover  on  the  diffracted  and  back- 
scattered  pressure.  These  results  are  illustrated  in  Chapter  III 
and  are  discussed  In  detail  In  Section  4.3. 

4. 2 Physical  Interpretation  of  the  Solution 

The  solution  for  the  diffraction  of  sound  radiated  by  various 
types  of  acoustic  sources  by  a locally  reacting  Impedance  covered 
half-plane  provides  an  opportunity  to  study  the  nature  of  the 
diffraction  phenomenon  as  well  as  the  effect  of  many  parameters. 
Since  there  are  three  types  of  acoustic  sources  considered  in  this 
study,  each  is  discussed  In  turn  emphasizing  the  differences  between 
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The  solution  for  Che  diffraction  of  a plane  wave  consists  of 
two  terms.  This  Is  evident  from  Che  solution  shown  In  Equation 
(2.88).  The  form  of  each  term  Is  really  very  similar,  differing 
slightly  In  three  places.  The  first  two  differences  are  the  plus  and 
minus  signs  in  front  of  In  both  the  exponents  and  In  Che  arguments 
of  the  Fresnel  Integrals.  The  third  difference  Is  In  Che  functions 

and  Each  of  the  two  terms  have  the  exponential  functions 

Ikrcos  (1i+6q)  -ik.rcos((})-e-)  hi  t u u i 

e ° and  e ° » respectively.  Taken  by  themselves, 

these  two  exponents  represent  plane  waves  propagating  In  the  direction 

of  Che  reflection  and  In  Che  direction  of  Incidence.  The  fact  that 

these  two  factors  are  present  In  Che  solution  hints  at  Che  existence 

of  both  a source  and  an  Image  source  In  Che  solution.  IncldenCly, 

Sommerfeld's  formulation  relies  heavily  on  the  concept  of  sources  and 

Image  sources.  These  two  terms  In  Che  diffraction  solution  are  not 

sources  and  Images,  however,  because  each  is  multiplied  by  a term 

containing  the  Fresnel  Integrals.  The  functions  containing  the 

Fresnel  Integrals  have  certain  Interesting  properties  that 

characterize  Che  diffraction  solution. 

When  Investigating  the  behavior  of  Che  Fresnel  Integrals,  It  Is 
most  Important  Co  consider  the  arguments.  The  arguments  of  the  two 
terms  are,  respectively: 
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and 

, 

►'  2kr  cos  — — . (^.2) 

1/2 

It  Is  Immediately  evident  that  each  argument  has  the  same  (kr) 

1/2 

dependence.  This  dependence  on  (kr)  manifests  Itself  In  two 

aspects  of  the  solution.  First,  as  kr  Increases  both  Fresnel 

Integrals,  C(X)  and  S(X)  approach  1/2  at  a rate  proportional  to 
1/2 

(kr)  . As  a result,  the  entire  quantity  that  contains  the  Fresnel 

-1/2 

Integrals  decays  to  zero  at  a rate  proportional  to  (kr) 

Subsequently,  the  diffracted  pressure  should  also  decay  by  the  same 

rate  which  Is  the  expected  decay  of  line  source  radiation.  Secondly, 

the  quantity  kr  Influences  the  oscillations  In  the  diffracted  pressure 

because  of  the  varying  phases  of  the  two  terms  of  the  solution. 

Because  of  the  nature  of  the  Fresnel  Integrals,  the  phase  varies  at 

a rate  proportional  to  the  square  of  the  argument.  This  means  that 

1/2 

the  phase  Is  directly  related  to  kr  rather  than  (kr)  . This  can 
be  demonstrated  by  comparing  Figures  3.1  and  3.2.  The  number  of 
oscillations  occurring  In  the  plots  Increases  In  direct  proportion 
to  kr.  As  kr  Is  doubled  from  kr  ■ 30  to  kr  ~ 100,  the  number  of 
oscillations  also  doubles. 

The  arguments  of  the  Fresnel  Integrals  In  Equations  (4.1)  and 
(4.2)  also  show  a dependence  upon  the  angle.  Note  that  at 
the  argument  In  Equation  (4.1)  goes  to  zero  while  the  other  argument 
given  In  Equation  (4.3)  goes  to  zero  at  (fi  ■ -iT+(p^,  It  Is  no  accident 
that  these  angles  correspond  to  the  reflection  showdow  boundary  and 
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the  incident  shadow  boundary,  respectively.  The  Fresnel  Integrals 

are  zero  when  the  arguments  are  zero.  As  a result,  the  two  terms 

that  contain  the  Fresnel  Integrals  are  maxima  at  their  respective 

shadow  boundaries.  This  behavior  of  the  arguments  contributes  to 

the  ability  of  the  diffraction  solution  to  accommodate  the  dramatic 

transition  that  occurs  across  the  shadow  boundaries.  Consider  for 

the  moment  only  the  first  term  in  Equation  (2.86)  evaluated  at 

ij)  ■ tt-A  . This  term  becomes 
o 


P (sin9^  - 
o 

2 (sin9^ 


sini))  ) e 
o 

+ sin())  ) 
o 


-ikr 


(4.3) 


Meanwhile,  the  argument  of  the  Fresnel  Integral  in  the  second  term 

is  / 2kr  sinij)^  and,  if  kr  is  sufficiently  large,  this  second  term 

-1/2 

as  a whole  will  be  negligible  because  of  the  (kr)  decay.  All 
that  survives,  then,  at  the  reflection  shadow  boundary  is  what 
appears  in  Equation  (4.3).  Incidentally,  this  value  has  the  specific 
physical  interpretation  of  being  reflected  wave  where 


(8ln<(!^  - sln9^) 
(sini})^  + 8in9^) 


(4.4) 


The  quantity  is  the  coefficient  of  reflections  of  a plane  wave 
incident  at  an  angle  on  an  Impedance  covered  infinite  plane  whose 
Brewster  angle  if  &■*■.  At  the  incident  shadow  boundary,  -Tr+())^,  the 
solution  equals  -p^/2.  This  Is  one-half  the  incident  wave.  The 
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behavior  at  the  two  shadow  boundaries  demonstrates  that  the  pressure 
Is  continuous  across  them.  At  the  dramatic  transition  between  the 
Illuminated  region  and  the  shadow  region,  the  pressure  becomes  the 
average  of  the  pressure  between  the  two  regions. 

It  has  been  previously  mentioned  that  the  figures  in  Chapter  III 
illustrate  dramatic  oscillations.  This  behavior  results  from  the 
Interference  of  waves  as  they  propagate  away  from  the  edge.  The 
exact  mechanism  that  causes  this  phenomena  Is  not  Immediately  clear 
from  the  solution.  To  explain  this  interference  phenomenon,  It  is 
helpful  to  consider  a plot  for  the  diffraction  of  an  incident  plane 
wave  at  ■ 120°  by  a perfectly  rigid  half-plane.  This  plot  is 
now  shown  in  Figure  4.1  for  two  different  values  of  kr.  Note  that 
there  exists  the  same  kind  of  oscillations  that  occur  for  an 
impedance  covered  half-plane.  This  suggests  that  the  oscillations 
do  not  result  from  the  influence  of  an  arbitrary  impedance  cover. 

If,  however,  the  absolute  values  of  the  two  terms  in  the  solution 
for  the  diffracted  pressure  are  plotted  separately  they  appear  as 
smooth  functions  as  illustrated  in  Figure  4.2. 

Note  that  these  two  smooth,  pear-shaped  plots  have  their 
respective  maxima  at  the  location  of  the  two  shadow  boundaries  and 
then  each  diminishes  away  from  that  angle  in  both  directions.  It 
seems  then  that  the  oscillations  results  from  the  interference 
between  these  two  terms.  Since  they  are  both  complex  functions,  the 
phase  differences  gives  rise  to  the  interference  behavior  as  they 
add  in  and  out  of  phase.  In  addition,  the  plots  in  Figure  4.2 


suggest  that  the  oscillations  result  from  the  wave  interference 
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developed  between  the  four  boundaries:  the  two  shadow  boundaries 
at  and  at  and  the  physical  boundaries  ± it.  In  light  of 

this  Interpretation,  It  Is  Interesting  to  look  at  a plot  for  each 
of  the  two  terms  when  the  angle  of  Incidence  matches  the  Brewster 
angle.  For  this  case,  there  should  be  no  reflection.  Figure  4.3 
represents  the  two  terms  for  a 120°  Incident  plane  wave  where  O"*"  • 60° 
and  6 “ 0°.  Note  that.  In  the  plot  for  the  first  term,  there  Is  a 

zero  at  the  exact  reflection  shadow  boundary.  However,  even  though 
there  Is  no  reflection,  there  Is  still  a contribution  away  from  the 
reflected  shadow  boundary.  The  plot  Is  almost  omnidirectional 
except  for  the  null  at  ^“0^-  On  the  other  hand,  the  plot  for  the 
second  term  clearly  Indicates  the  shadow  boundary  at  -11+4)^.  Again 
the  interference  phenomenon  occurs  when  these  two  complex  terms 
are  summed  and  plotted  In  Figure  4.3.  In  short,  the  Implication 
Is  that  the  oscillations  In  the  dlf'  acted  pressure  are  a result 
of  the  Interference  between  the  waves  Interacting  with  the  four 
boundaries,  the  two  shadow  boundaries,  and  the  two  physical  boundaries 
at  t TT. 

The  fact  that  at  the  exact  shadow  boundary  the  diffracted 
pressure  converges  to  the  mean  value  between  the  two  regions  Is 
particularly  Interesting  when  viewed  In  comparison  with  the  result 
given  by  the  far-fleld  solution  at  the  exact  shadow  boundary.  The 
far-fleld  solution  Is  given  In  Equation  (2.49).  In  the  far-fleld 
solution  there  are  singularities  at  4 “ at  <})  " 

because  of  the  denominator.  Consequently,  the  solution  blows  up 
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at  the  shadow  boundaries.  This,  of  course,  Is  not  realistic.  One 
does  not  actually  have  an  infinite  pressure.  To  explain  the  apparent 
fallacy,  recall  that  the  far-field  solution  results  from  the  steepest 
descent  approximation.  If  the  saddle  point,  however,  lies  near  a 
singularity,  the  method  falls.  This  is.  Indeed,  the  case  when 
$ “ or  i}>  “ Consequently,  the  far-field  solution  is  not 

really  valid  when  it  predicts  infinite  pressure  at  the  shadow 
boundaries.  It  is  in  these  regions,  close  to  the  shadow  boundaries, 
where  the  behavior  of  the  Fresnel  Integrals  is  necessary  to  account 
for  the  dramatic  transition  across  the  shadow  boundaries. 

It  is  appropriate,  though,  in  cases  where  kr  is  very  large  to 

accept  the  far-field  solution  because  it  does  Illustrate  that  the 

greatest  diffraction  effect  does  occur  at  the  shadow  boundaries. 

It  is  good  enough  to  satisfy  one's  physical  intuition  so  long  as 

-1/  2 

one  recalls  that  the  pressure  has  diminished  at  a rate  of  (kr) 
and  that,  at  the  exact  boundary,  what  appears  to  be  an  Infinity 
should  really  be  one-half  the  mean  between  the  two  regions. 

The  far-field  plots  are  helpful  in  searching  for  the  nulls  in 
the  diffracted  and  backscattered  fields.  If  the  solution  is 
considered  in  the  form  given  in  Equation  (2.49),  an  expression  for 
the  zeroes  of  the  diffracted  pressure  and  backscattered  pressure 
may  be  found.  These  expressions  may  be  found  in  Equations  (2.31) 
and  (2.32). 

Much  of  what  occurs  in  the  plane  wave  solution  also  occurs  in 
the  solutions  for  the  line  source  and  the  point  source.  All  three 
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are  essentially  examples  of  constructive  and  destructive  wave 
interference  between  the  two  shadow  boundaries  and  the  two  physical 
boundaries.  There  are,  however,  a number  of  minor  differences  in 
each  solution  that  require  some  amplification. 

The  dlf traction  solution  for  a line  source  is  given  in  Equation 
(2.113).  One  can  apply  a limited  geometric  Interpretation  to  this 
solution  by  once  again  thinking  in  terms  of  a line  source  and  a 
corresponding  image.  In  order  to  appreciate  this,  consider  Figure 
4.4. 

The  parameters  that  appear  in  the  solution  are  R^,  R,  and  S. 

From  Figure  4.4,  the  values  R and  S are  the  distances  from  the 
source  and  image  to  the  observer  point,  while  the  value  R^^  « r + r^ 
is  the  path  distance  a wave  must  travel  from  the  source  to  the  edge 
and  from  there  to  the  observer. 

The  famous  reciprocity  theorem  is  also  verified  by  this  solution. 

If  the  values  for  kr  and  kr^  are  reversed  in  Equation  (2.113),  which 

is  essentially  reversing  the  position  of  the  source  and  observer, 

there  is  no  change  in  the  solution. 

The  solution  for  the  point  source  is  given  in  Equation  (2.132). 

Like  the  solution  for  the  line  source,  there  is  a similar  geometric 

interpretation.  Also  note  that  for  this  solution,  the  reciprocity 

theorem  is  once  again  validated.  By  reversing  the  values  of  kr 

and  kr  and  also  the  values  for  kz  and  kz  , the  result  remains 
o o 

unaltered. 

In  the  point  source  solution,  Che  variable  kR^  has  an 
interesting  physical  interpretation.  The  length  R^  is  Che  shortest 

» 

i 
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distance  a ray  may  travel  from  the  source  to  the  edge  of  the  half- 
plane and  finally  to  the  observer.  The  appearance  of  the  shortest 
path  distance  in  the  formula  Is  a manifestation  of  Fermat's  principle 
on  the  dlftraction  of  a bundle  of  rays.  The  fact  that  appears  in 
the  solution  verifies  the  principle.  Figure  4.5  Illustrates  Che 
geometry  for  the  point  source,  image  point  source  and  observer.  The 
dashed  line  indicates  Che  shortest  distance  a ray  would  travel  when 
going  from  the  source  to  the  edge  and  finally  to  the  observer.  The 
geometry  for  the  path  of  the  ray  is  detailed  in  Figure  4.6.  From 
this  figure,  one  realizes  that  the  triangles  OAB  and  SAC  are  similar, 
the  length  from  B to  C is  (z-z^),  the  length  OB  is  r and  the  length 

SC  is  r . It  follows  that 
o 

(4.5) 

r r 

o 


and 


(z-z  ) • AB  + AC 
o 

r 


- AB  + AB  — 

r 

r+r 


Consequently 


AB 


r(z-z  ) 


(4.6) 


(4.7) 
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Figure  4.5. 


The  geometry  for  the  source,  image  source,  and 
observer  for  an  Incident  line  source. 
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and 


AC 


r (z-z  ) 
o o 

(r+r  ) 


(4,8) 


The  distance 


, z-z  2 

^ o 

i hr 


and  the  distance 


(4.9) 


AS 


z-z 

o 


r+r 

o 


2 


(4,10) 


Finally,  one  can  calculate  the  distance  from  Equations  (4.9) 
and  (4.10)  as  follows: 


= OA  + AS 


r^I(r+r  )^  + (z-z  )^]  r ^[(r+r  )^  + (z-z  )^ 

o o ^ _o o o 


(r+r  )‘ 
o 


(r+r  )‘ 
o 


/Tr+r  )^  + (z-z 


(4.11) 


in  summary,  the  solution  for  the  diffraction  of  tlie  radiation 
of  each  type  of  acoustic  source  illustrates  the  occurrence  of  con- 
structive and  destructive  Interference  and  demonstrates  the  validity 
of  the  Keciproclty  Theorem  and  Fermat's  Principle. 
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4 . 3 Discussion  of  Results 

The  numerical  results  that  were  compiled  In  Chapter  111  are 
discussed  In  this  section. 

For  the  plane  wave  one  may  compare  the  transition  from  the 
near-fleld  solution  to  the  far-fleld  solution  from  any  of  the 
Figures  3.1  to  3.14.  For  example,  in  Figures  3.1  and  3.2,  it  is 
evident  that  increasing  kr  results  in  a more  complicated  phase 
Interaction  as  Indicated  by  the  Increase  In  the  number  of  oscilla- 
tions In  the  plots.  It  is  also  curious  that  even  at  kr  • 100, 
which  represents  approximately  15  wavelengths  from  the  edge,  the 
solution  does  not  yet  produce  the  clearly  defined  null  that  appears 
In  the  plot  of  the  far-fleld.  This  Is  a consequence  of  the  fact 
that  the  exact  phase  cancellation  required  to  form  this  null  only 
occurs  at  kr  equal  to  infinity,  which  does  not  exist  in  reality. 

The  first  eight  figures  show  how  the  diffracted  pressure 
varies  as  a function  of  the  impedance  cover  on  the  Illuminated 
surface.  Specifically,  compare  the  Figures  3.2,  3.4,  3.6,  and 
3.8  where  the  Brewster  angles  for  the  Illuminated  surface  are 
0+  - 30°,  60°,  90°,  and  i“  (pressure  release),  respectively.  It 
appears  that  as  this  Impedance  cover  on  the  Illuminated  surface 
varies,  there  is  a significant  change  in  the  diffracted  pressure. 

In  contrast,  compare  these  same  four  figures  with  Figures  3.10 
and  3.12,  where  the  Brewster  angle  on  the  unlllumlnated  surface  is 
varied  from  the  extremes  of  ■ 90°,  for  a totally  absorbing 
surface  to  ■ 51,  for  a pressure  release  surface.  What  becomes 
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evident  from  this  comparison  is  that  the  impedance  of  the  surface 
lying  in  the  shadow  region  has  little  effect  on  the  diffracted 
pressure  while  tne  Impedance  cover  on  the  illuminated  surface  has 
a significant  effect,  being  most  apparent  in  the  illuminated  region. 

One  concludes  that  the  Impedance  cover  of  the  illuminated  surface 
dominates  the  diffracted  pressure  in  the  illuminated  region.  Although 
the  surface  condition  on  the  opposite  side  may  Influence  the  diffracted 
pressure,  it  is  to  a much  smaller  degree. 

Figures  3.13  and  3.14  are  interesting  in  contrast  to  Figures 
3.1  and  3.2.  The  impedance  cover  of  the  half-plane  in  Figures  3.13 
and  3.14  is  almost  equivalent  to  that  in  Figures  3.1  and  3.2  except 
for  a small  imaginary  (reactive)  component  added  to  the  Impedance  on 
the  upper  surface.  This  small  imaginary  part  seems  to  have  little 
effect  on  the  near-field  solution.  However,  in  the  far-fleld,  where 
a true  null  is  shown  in  Figure  3.2,  only  a minimum  is  shown  at  approx- 
imately the  same  angle  in  Figure  3.14.  It  is  evident,  then,  that 
nulls  exist  in  the  far-field  solution  only  when  the  impedance  is 
real,  i.e.,  absorptive.  In  addition,  if  the  Impedance  is  complex, 
with  a small  imaginary  part,  the  null  becomes  a minimum  at  approx- 
imately the  same  location.  Also,  increasing  the  imaginary  part 
further  may  remove  the  minimum  altogether. 

The  plots  of  the  backscattered  pressure  further  establish  the 


conclusion  that  the  effect  of  the  Impedance  on  the  Illuminated 
surface  dominates  the  backscattered  pressure  in  the  half-space  that 
faces  the  illuminated  surface.  Compare  Figures  3.16,  3.18,  3.20, 
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and  3.22.  In  each,  the  lower  surface  Is  rigid  while  the  Brewster 
angles  for  the  upper  surface  are  O"*”  = 30°,  60°,  90°  and  15,  respec- 
tively. Recall  that  for  backscattered  pressure,  the  source  and 
observer  are  at  the  same  location.  Then,  the  backscattered  pressure 
In  the  lower  half-plane  appears  to  be  the  same  for  all  of  the 
Brewster  angles  0 since  the  illuminated  lower  surface  Is  always 
rigid.  The  backscattered  pressure  remains  the  same  In  the  lower 
half-plane  regardless  of  what  happens  on  the  upper  surface.  In 
contrast,  however,  when  the  source  and  the  receiver  are  in  the  upper 
half-plane  the  backscattered  pressure  from  the  different  plots 
changes  corresponding  to  each  different  Impedance  conditions.  The 
Implication  Is  that  the  Impedance  condition  on  the  llluialnated 
surface  has  the  dominant  influence  on  the  behavior  of  the  back- 
scattered  pressure. 

A comparison  of  the  plots  for  both  the  line  source  and  the 
point  source  to  those  for  the  plane  wave  incidence  demonstrates  a 
similar  dependence  on  the  Impedance  cover. 

4 . 4 Summary  and  Conclusions 

The  problem  of  diffraction  and  backscattering  of  various  types 
of  sources  was  attempted  by  means  of  a so-called  dual  Integral 
formulation.  The  dual  Integral  formulation  resulted  in  a new  and 
closed  form  solution  for  the  diffraction  of  various  acoustic  sources 
by  an  Impedance  covered  half-plane.  The  closed  form  solution  illus- 
trates that  the  diffraction  is  a phenomenon  manifested  by  wave 


Interference. 
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Of  an  exhaustive  number  of  numerical  results  compiled  for  many 
parameter  combinations,  only  a portion  are  presented  In  this  thesis. 
They  demonstrate  that  the  effect  of  the  surface  Impedance  dominates 
in  the  Illuminated  half-space  Into  which  it  faces,  even  though  the 
influence  of  each  surface  extends  with  diminishing  effect  Into  the 
opposite  half-space  behind  the  barrier. 

At  the  geometric  shadow  boundary,  the  diffracted  pressure  becomes 
one-half  the  Incident  pressure,  while  at  the  reflection  shadow 
boundary,  the  diffracted  pressure  Is  equal  to  one-half  the  reflected 
pressure  from  an  Impedance  covered  Infinite  surface. 

The  results  also  Indicate  that,  although  the  near-fleld  solution 
approaches  the  far-field  solution  as  kr  becomes  large.  It  never 
equals  the  far-fleld  solution.  Both  the  source  and  observer  must  be 
at  an  infinite  distance  from  the  edge  to  allow  the  precise  phase 
cancellations  necessary  to  yield  the  far-fleld  solution. 
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